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Abstract 

We study the thermodynamics of maximally supersymmetric U (N) Yang-Mills theory 
on R X S*^ at large A^. The model arises as a consistent truncation of A/" = 4 super 
Yang-Mills on R x S''^ and as the continuum limit of the plane-wave matrix model ex- 
panded around the spherical membrane vacuum. The theory has an infinite number 
of classical BPS vacua, labeled by a set of monopole numbers, described by dual super- 
gravity solutions. We first derive the Lagrangian and its supersymmetry transformations 
as a deformation of the usual dimensional reduction of A/" = 1 gauge theory in ten di- 
mensions. Then we compute the partition function in the zero 't Hooft coupling limit in 
different monopole backgrounds and with chemical potentials for the i?-charges. In the 
trivial vacuum we observe a first-order Hagedorn transition separating a phase in which 
the Polyakov loop has vanishing expectation value from a regime in which this order pa- 
rameter is non-zero, in analogy with the four- dimensional case. The picture changes in 
the monopole vacua due to the structure of the fermionic effective action. Depending on 
the regularization procedure used in the path integral, we obtain two completely differ- 
ent behaviors, triggered by the absence or the appearance of a Chern-Simons term. In 
the first case we still observe a first-order phase transition, with Hagedorn temperature 
depending on the monopole charges. In the latter the large behavior is obtained by 
solving a unitary multi-matrix model with a peculiar logarithmic potential, the system 
does not present a phase transition and it always appears in a "deconfined" phase. 
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1 Introduction 



In the context of the AdS/CFT correspondence [H, [2], [3], H] an interconnected family 
of theories with sixteen supercharges has been recently studied [5]. They all have a 
mass gap and a discrete spectrum of excitations. These theories can be obtained from 
consistent truncations of A/" = 4 super Yang-Mills on R x 5*^ and have many BPS vacua. 
Remarkably, smooth gravity solutions corresponding to all these vacua can be described 
rather explicitly. At large 't Hooft coupling some properties of the dual string theory have 
also been examined according to the pioneering proposal of p]. 

From the gauge theoretical point of view it seems particulary appealing to investigate 
the properties of one specific theory belonging to this class, namely the maximally su- 
persymmetric U{N) Yang-Mills theory on R x S"^. This theory already appeared in [7] 
where it arises from the fuzzy sphere vacuum (membrane vacuum) of the plane-wave ma- 
trix model by taking a large limit that removes the fuzzyness. The model can also be 
constructed from the familiar A/" = 4 SYM theory by truncating the free-field spectrum 
on R X 5^ to states that are invariant under U{1)l C SU{2)l, where SU{2)i is one of the 
SU{2) factors in the 5*0(4) rotation group of the three-sphere. Geometrically this corre- 
sponds to a dimensional reduction of the four-dimensional supersymmetric theory along 
the U{1) fiber of S'^ seen as an Hopf fibration over S"^. The resulting model lives in one di- 
mension less and maintains supersymmetry through a rather interesting mechanism. The 
particular dimensional reduction breaks the natural 5*0(7) i?-charge symmetry to 5*0(6), 
singling out one of the seven scalars of the maximally supersymmetric Yang-Mills theory, 
which then behaves differently from the others. It combines with the gauge fields to form 
a peculiar Chern-Simons-like term that is crucial to preserve the sixteen supercharges, 
balancing the appearance of mass terms for fermions and scalars. The BPS vacua are 
generated by the same term that allows to combine the field strength and the scalar into 
a perfect square whose zero-energy configurations are determined by A^ integers ni, tt-tv 
associated to monopole numbers on the sphere. 

The model represents an interesting example of a supersymmetric non-conformal gauge 
theory, with smooth gravitational dual and non-trivial vacuum structure, defined on a 
compact space. The last feature is particulary appealing in the study of the thermal 
properties of the theory. Recently the thermodynamics of large A^ theories on compact 
spaces has attracted much attention. On compact spaces the Gauss's law restricts physical 
states to gauge singlets. Consequently, even at weak 't Hooft coupling the theories are 
in a confining phase at low temperature and undergo a deconfinement transition at a 
critical temperature. For example, the partition function of A/" = 4 super Yang-Mills 
theory on R x 5^ was computed at large A^ and small coupling in [HI [9], [10]. It was 
shown that the free energy is of order 0{1) at low temperature and of order 0{N'^) 
above a critical temperature. At strictly zero 't Hooft coupling the transition is a first- 
order Hagedorn-like transition. At small coupling a first or a second order transition is 
expected, depending on the particular matter content of the theory. The computation in 
the Af = A maximally supersymmetric case has never been performed but in [11] it was 
argued that the maximally supersymmetric plane- wave deformation of Matrix theory and 
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A/" = 4 SYM should show similar behavior, including thermodynamics. The plane wave 
matrix model is a theory with sixteen supercharges and it was argued in |^ to be dual 
to a little string theory compactified on S^. For a small sphere, this theory is weakly 
coupled and one may study the little string theory thermodynamics rather explicitly |12j . 
The phase transition for this model was shown to remain first order in [13] indicating that 
this might also be the case for Af = A SYM. This was shown by computing the relevant 
parts of the effective potential for the Polyakov loop operator to three loop order |13j . 
With the same procedure it was shown in [T3] that also for pure Yang-Mills the phase 
transition remains first-order up to three loops. The phase transition at weak coupling 
is basically driven by a Hagedorn-like behavior of the spectrum in the confining phase, 
suggesting a possible relationship with the dual description of large N gauge theories 
in terms of strings. For A/" = 4 the relevant string theory lives on an asymptotic AdS 
space and, at large 't Hooft coupling, the deconfinement phase transition corresponds to a 
Hawking-Page transition [151 IS] . The thermal AdS space dominates at low temperature 
and the AdS-Schwarzschild black hole is the relevant saddle-point in the high-temperature 
regime. The original proposal presented in [SI (H] to connect the phase transitions at small 
coupling on compact spaces with the gravitational/stringy physics stimulated a large 
number of investigations. Lower- dimensional theories on tori were examined in [TTJ [18], 
while the inclusion of chemical potentials for the i?-charges was discussed in [191 120] and, 
more recently, pure Yang-Mills theory on [2T] was found to have a second order phase 
transition at small 't Hooft coupling. 

In this paper we study the thermodynamics of A/" = 8 super Yang-Mills theory on 
]R X S*^. We first derive the Lagrangian and its supersymmetry transformations as a de- 
formation of the usual dimensional reduction of A/" = 1 gauge theory in ten dimensions. 
Actually our procedure will generate a larger class of three-dimensional theory: according 
to the particular choice of the generalized Killing spinor equation we obtain also theories 
on AdSa with peculiar Chern-Simons couplings. Then we compute the A/" = 8 parti- 
tion function in the zero 't Hooft coupling limit, for different monopole vacua. In the 
trivial vacuum we observe a first-order Hagedorn transition separating a phase in which 
the Polyakov loop has vanishing expectation value from a regime in which this order pa- 
rameter is non-zero, in complete analogy with the four- dimensional case. The Hagedorn 
temperature is also obtained in the presence of chemical potentials for the -R-charges. 
Discussions on the dual gravitational picture [5] and the possibility of matching the gauge 
theory Hagedorn transition with a stringy Hagedorn transition, by exploiting for example 
a decoupling limit as in [201 1221 [231 121] postponed to a forthcoming investigation. 

The situation is very different in the non-trivial monopole vacua. The original U (N) 
gauge group is broken to a direct product U{Ni) x U{N2) x ..U{Nk) and the constituent 
fields transform, in general, under bifundamental representations of U{Ni) x U{Nj). Be- 
cause of the Gauss's law on a compact manifold, however, the only allowed excitations 
are SU{Ni) x SU{Nj) singlets. Different selection rules are instead possible for the U{1) 
charges in three dimensions, depending on the definition of the fermionic Fock vacuum 
in the presence of background monopoles [25] . The appearance of fermionic zero-modes 
makes possible, in general, to assign a non-trivial charge to the Fock vacuum, as clearly 
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explained in [26]. In the path- integral formalism this corresponds to precise choices in 

regularizing fermionic functional determinants which might produce Chern-Simons terms 
in the effective action. In our case the different possibilities are clearly manifested in the 
matrix model describing the partition function. We recall that, in the trivial vacuum, the 
thermal partition function is reduced to an integral over a single U (N) matrix [HI [S] 



where U = e*'^" (a is the zero mode of the gauge field Aq on x and (3 = \/T the 
inverse of the temperature). In the non-trivial monopole vacuum Z{I3) is given instead 
by a multi-matrix model over a set of unitary matrices Ui{Ni), i = 1,2, ..k, reflecting the 
breaking of the U{N) gauge group. More importantly the effective action SeffiUi), at zero 
't Hooft coupling, can be modified by the presence of logarithmic terms NQiTi \og{Ui) 
that implement selection rules on the U{1) charges. The large analysis is highly affected 
by these new interactions: they contribute at order A^^ and can drive the relevant saddle- 
point always at a non-zero value of the Polyakov loop. Unitary matrix model of the kind 
we encountered in our analysis have been previously considered in the eighties [271 [25] , but 
with an important difference: in those studies the coefficient weighting the logarithmic 
term Tr \og{U) in the action was taken independent on A^. Conversely the large N saddle- 
points were not modified by its presence, being determined by the rest of the action. In 
our case, instead, we have to cope with a linear dependence on and we cannot simply 
borrow those results. We have therefore performed an entirely new large A^ analysis of 
these kind of models, starting from an exact differential equation of the Painleve type 
that describes the finite A^ partition function [29] . 

The paper is organized as follows. In section 2 we construct the supersymmetric 
Yang-Mills theory on R, x S*^ using a different strategy with respect to [5| and [7| (see 
also [3U| for a careful derivation of the Hopf reduction and [HT] for an extension to more 
general fiber bundles). We start from M = 1 super Yang-Mills theory in ten dimensions 
and consider its dimensional reduction on IR x 5*^. We find the relevant Killing spinors 
that generate the rigid supersymmetry, generalizing to our case the approach developed 
in [32| . We further determine the deformations of the original ten dimensional Lagrangian 
and of the supersymmetry transformations ensuring the global invariance of the action. 
Interestingly, using the same strategy it is possible to construct two other maximally 
supersymmetric gauge theories on three-dimensional curved spacetimes, living both on 
AdSs and differing from the theory introduced in [7| in the structure of the Chern-Simons 
terms. In section 3 we briefly examine the BPS vacua of the model, we comment on their 
gravitational description and the related instanton solutions. 

We then turn to study the thermodynamics at zero 't Hooft coupling. Following the 
analysis in [HI [H|, we obtain the partition function of the theory in a generic vacuum, in 
terms of matrix integrals. In section 4 we present the results of the relevant functional 
determinants in the background of a gauge flat- connect ion and of a monopole poten- 
tial, recovering the appropriate single-particle partition functions for scalars, spinors and 
vectors. Careful ^-function evaluations are deferred to the appendices. We discuss the 
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emerging, on the monopole background, of new logarithmic terms in the effective action, 
directly related, in this formalism, to the appearance of fermionic zero-modes. We explain 
their dependence on the regularization procedure and remark their interplay with a typ- 
ical three-dimensional phenomenon, the induction of Chern-Simons terms. We interpret 
their effect as a part of the projection into singlets of the gauge group, as required by the 
Gauss's law. Section 5 is devoted to discuss the large thermodynamics in the trivial 
vacuum. We determine the critical temperature at which the first-order phase transition 
takes place and we generalize the result to the case of non vanishing chemical potentials 
for the i?-charges. Finally, in sections E] and [TJ we study the large theory on the non- 
trivial monopole backgrounds: we consider a large class of vacua, characterized by the 
set of integers ni,..,nk and large degeneracies Ni,..,Nk. According to the discussion 
of section 4, we study two different choices for the logarithmic terms, within our regu- 
larization procedure. First, in section [6l we discuss the "uncharged" case, that amounts 
to make a particular choice of branch cuts, in the ^-function regularization procedure 
[331 El], that cancels the Chern-Simons like contributions. In turn we get a non- vanishing 
Casimir energy, depending explicitly on the monopole background. The resulting uni- 
tary multi-matrix model is an obvious generalization of the trivial case. We find again a 
first-order phase transition, with an Hagedorn temperature explicitly depending on the 
monopole numbers. We discuss also some particular class of vacua, characterized by large 
monopole charges, whose Hagedorn temperature approaches the one of the theory on 
S^/Zk in trivial vacuum. In section [7] we discuss the opposite situation of a "maximally" 
charged fermionic vacuum: we have a non-trivial modification of the unitary multi-matrix 
model due to appearance of the new logarithmic terms and vanishing Casimir energy. 
For the sake of clarity we will restrict our discussion to a particular simple background 
[n, n, ..,n, —n, —n.., —n). We show the existence of a non-trivial saddle-point for the ef- 
fective action for a wide range of temperatures starting from zero, within the assumption 
that we can disregard higher windings contributions in this regime. This implies that the 
theory is always in a "deconfined" phase. We have to face the problem of computing the 
free energy and the phase structure of the matrix model 



that is a non-trivial deformation of the familiar Gross- Witten model [3S|- Its large A^ 
behavior is carefully studied in section 7.1 , obtaining the exact free energy in terms of 
the solution of a fourth-order algebraic equation: we prove that there is no phase transition 
as long as p 7^ 0, in contrast with the usual p = case, that appears as a singular point 
in the parameter space. In section 7.2 we use the results of our analysis to derive a set of 
saddle-point equations for the partition function which describes the "deconfined" phase. 
The disappearance of the confining regime is consistent with the known results on finite 
temperature 2+1 dimensional gauge theories where, once a topological mass (a Chern- 
Simons term) is turned on, there cannot be a phase transition [361 EZl EH]. In section 8 
we briefly draw our conclusions and discuss future directions. Several appendices are 
devoted to technical aspects and to an alternative derivation of the partition functions. 
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In appendix A we report some details on supersymmetry transformations. In appendix 
B we give the details of the computation of functional determinants. In appendix C 
we recover the results for the single-particle partition functions from those of the parent 
A/" = 4 theory by explicitly constructing the projector into the U{1) invariant modes. We 
also check the consistency of our results with those of [SH], where the theory on R x S'^/Z^ 
has been studied. Appendix D is instead focused on some technical aspects, related to 
the solution of the large matrix integrals. 

2 Lagrangian and supersymmetry on R x 5^ from 

D = 10 

There are many ways to construct the Lagrangian of the gauge theory with sixteen super- 
charges on R X S*^ and its supersymmetry transformations. For instance, in [7] this theory 
was obtained from the plane- wave matrix model action expanded around the /c-membrane 
vacuum in the large A^ limit. Subsequently, in [5] it was derived as a U{1) truncation 
of the spectrum of the A/" = 4 gauge theory on IR x S"^. Since here we shall be mainly 
concerned with the field theoretical features of this A/" = 8 model, we shall follow a more 
conventional (and maybe pedagogical) approach: the Lagrangian and its supersymmetry 
transformations will be derived as a deformation of the standard toroidal compactification 
of A/" = 1 gauge theory in ten dimensions. 

We first consider the theory on the flat Minkowski space in three dimensions, M(i 2). 
The A/" = 8 theory in this case is the straightforward dimensional reduction of the A/" = 1 
theory in D = 10. The most convenient and compact way to present its Lagrangian is 
to maintain the ten-dimensional notation and to write (see appendix |X] for a summary of 
our convention^ 

C''' = + ti^T^'DMi^. (2.1) 

All the fields in (12.11) only depend on the space-time coordinates (a;°, x^, x^). In particular, 
from the three-dimensional point of view, the gauge field Am contains the reduced gauge 
field and seven scalars (0^) = (03, 04, ■ ■ ■ , ^9) = (03, 0m)- The flat ten dimensional 
space-time metric is diagonal and it has the factorized structure T'' x lM(i 2). 

Our goal is now to promote the supersymmetric theory in the flat 2+1-dimensional 
space-time to a supersymmetric theory on the curved space IR x 5*^. It is useful to 
keep a ten-dimensional notation where the above space-time is viewed as a submanifold 
embedded in x R x 5*^ with the metric 

7 

ds'' = -de + R'{de'' + sin' ed^^) + dr)^ . (2.2) 

1=1 

Here the coordinates 6 and span the sphere S*' of radius R, while the internal angular 
coordinates rji parameterize the torus T''. The action (12. ip in the background (12. 2p is 

^In general we shall omit the trace over the gauge generator in our equations, unless it is source of 
confusion. 
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still meaningful once we introduce the appropriate dependence on the vielbein and the 
spin-connections in the covariant derivatives. The real issue is whether this theory will 
have any supersymmetry. The action (12.11) on flat space is invariant under the usual 
supersymmetry transformations written in terms of a constant arbitrary spinor e 

6^°^ Am = -2iiljrMe, 

Constant spinors however do not exist, in general, on a curved space. For a space-time of 
the type fl2.2l) . the notion of a constant spinor should be replaced with that of a Killing 
spinor [32] • Its specific definition may depend on the detail of the geometry, but, for us, 
it will be a spinor satisfying an equation of the type 

V^e = K;T,T'^\ , (2.4) 

where the Greek indices run only over the three-dimensional space-time since the trans- 
verse coordinates rji are fiat and we can always choose e to be a constant along these 
directions. In (12.41) we have also inserted an additional dependence on the F matrices 
through a monomial factor F^^^ [§. This has double role: (a) it makes (12. 4p compatible 
with the ten-dimensional chirality conditions; (b) it generates, as we shall see, the rele- 
vant massive deformations for our fields. Finally the tensor ^ expresses an additional 
freedom in constructing the Killing spinors. In a curved space, there is in fact no a priori 
reason to treat all the coordinates symmetrically. In the IR x S"^ curved space-time ge- 
ometry there is a natural splitting between space and time and thus it is quite natural to 
weight them differently by choosing 

K^'^ = a[{6; + k,k^)-Bk,k'^], (2.5) 

where k^ is the time-like Killing vector of (12. 2p and a, B are two arbitrary parameters. 
The parameter a is fixed by imposing the necessary integrability condition (the first) |40j . 
which arises from the commutator [V^, V^je. This can be either expressed in terms of the 
space-time curvature scalar 7?. = 2/i?^ or, through (12. 4p . in terms of and consequently 
of a. We thus get for a 

The parameter B, instead, remains free and it will be determined in the following. 

The variation of the action (12. ip with respect to the supersymmetry transformations 
(12. 3p written in terms of a non-constant supersymmetry parameter e does not vanish. 
Terms depending on the covariant derivatives of e (12. 4p are in fact generated (see appendix 



^The direction (1, 2) span the tangent space to the sphere S'^, while the index 3 is along the first of 
the compactified dimensions. 
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A for conventions and more details) 
S'''C''' =2Re{ti^FMNm^'''V,e} 

(2.7) 

where in the second equahty we have used (12.41) and (12.51) . This undesired variation can 
be compensated by adding the following deformations to the original Lagrangian 

= iMa^T^^^ip + Na(f)3Fi2, = Va'^(j)'i + Wa^(f)l (2.8) 

and by adding new terms to the supersymmetry transformations of the fermions 

5''V = P«r"^r^23^^^ ^ GaT'T'''<j)se, (2.9) 

where M, N, V, W, P, G are arbitrary parameters to be fixed by imposing the invariance 
of the complete action. The size of the deformations is tuned by the natural mass scale 
a = 1/ {2R) provided by the radius of the sphere. 

Some comments on the form of (12. 8p and (12.91) are in order. The addition of mass 
terms for the scalars {C ) is a common and well-known property for supersymmetric 
theories in a background admitting Killing spinors. Some of the mass terms can also 
be justified with the requirement that the conformal invariance originally present in flat 
space is preserved. In four dimensions, for A/" = 4 super Yang-Mills, this is the only 
required modification of the Lagrangian because of an accidental cancellation. Since we 
are in three dimensions, we are also forced to introduce a non-standard mass term for the 
fermions (the first term in £^^*). The natural supersymmetric companion for a fermionic 
mass in D = 3 is then a Chern-Simons-like term (the second term in £ ^ ). Its unusual 
form, 03-Fi2, mixes the scalar (j)^ with the gauge-fields and is inherited from the particular 
choice of the monomial F^^^ in (12.41) . Then the modifications (12.91) in the supersymmetry 
transformations are the only possible ones with the right dimensions and compatible with 
the symmetries of the theory. 

The most convenient and simple way to analyze the effect of the additional terms in 
the Lagrangian (12.81) and in the supersymmetry transformations (12.91) is to single out, in 
the variation of the Lagrangian, different powers of the deformation parameter a. We 
start with the linear order in a, the zeroth order being automatically absent since our 
theory is supersymmetric in fiat space-time. At this order we have three contributions: 
the original variation (12. 7p . the variation of the new Lagrangian C^^^ with respect to the 
old transformations (12. 3p 

5V=2MaRe{ii^{Fijr^- 2FoiT^'- 2Do03r°=^+2A03r^=^ + 2Do0mrO™- 2A</'7nr'^+ 
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and finally the variation of C^°^ with respect to (12.91) 

(1) (0) 

5 C =2-Re{ia^{PV^^D^(t)^-^gPV^^[<P^,(t)n\\GV^^D^^^-^gGV^^[<^^^^ (2.11) 

See appendix A for all the different index conventions. It is quite straightforward to derive 
(12.101) and (12.111) since at this order in a we can consider e as a constant spinor, namely 

(0) (0) (0) (1) (1) (0) 

V^e = 0. Imposing that 5 C +5 C +5 C = C (a ) gives a linear system of eight 
equations in the five unknowns M, N ,P ,G and B. The details are given in appendix 
\KJ[ Quite surprisingly, this system is still solvable and it fixes the value of the above 
constants as 

M = ~, N = 4:,P = -2, G = -2, ^ = \- (2-12) 

The next and final step is to consider the order in our supersymmetry variation. The 
situation is much simpler now since we need to evaluate only few terms. We have in fact 
to consider the effects of the corrected transformation (12. 9p on C^^^ 

= iMad^'^ (^r^23^) ^ 2Re{za2^(r^0^ - 2T^<j)^)e) (2.13) 

and we have to take care of the terms coming from 5''^^ C''°^ originated from the covariant 
derivative of the Killing spinor e. We obtain 

d'^'C^"' = -2Re{za'^[3r™</.wr + QT^hV). (2.14) 

These two contributions are easily compensated by the variation of Z^'^' , 

5(o)£(2) = -Aia'{V(t>rni'T^i^ + {V + l^)03?r=^V^) , (2.15) 

By setting V = —1 and W = —3 no surviving term is left! We remark that there is 
no 0{a^) term, because there is neither an a-dependent term in the variation of bosons 
(which might produce a 0{a^) term in the variation of £ ^ ) nor term in the variation 
of fermions. 

We have thus reached our original goal: to promote the Af = 8 theory in flat space in 
three dimensions to an A/" = 8 theory in the curved background R x S*^. Its Lagrangian 
in a ten-dimensional language is thus given by 

C = -IfmnF'''' + iWDm^ - t^Wi^ + 2/.</>3Fi2 - ^01. - /iVL (2-16) 
and it is invariant under the supersymmetry transformations 
5Am = -2i'il)TMe, 

where fi is the mass-scale fi = 1/R. Notice that the mass for the scalars (with m = 
4, 5, . . . , 9) in (12.161) is that required by conformal invariance on R x S'^: ml^f ~ ~ 



= The mass of the scalar 03 is, instead, different because 03 mixes with the gauge 
fields. This mixing also breaks the original 50(7) i?-symmetry present in flat space 
to the smaller group 50(6)^ (~ 5'f/(4)/j): the bosonic symmetries IR x 5*0(3) x S0{6)r 
combine with the supersymmetries into the supergroup SU (2|4). We have to mention that 
our presentation heavily relies on the general analysis of [32] , where the problem of the 
existence of globally supersymmetric Yang-Mills theory on a curved space was addressed 
and some general recipes on how to construct these models were given. However, the 
Lagrangian (12.161) does not directly belong to the famihes of theories discussed in [32] , it 
realizes nevertheless a straightforward generalization of them. We have in fact allowed 
for a more general Killing spinor equation both by including the additional matrix factor 
K^"^ and by considering a monomial factor F^^^ mixing one of the transverse compact 
directions with the two spatial directions of the actual space-time of the theory. 
The Lagrangian fl2.16p written in terms of the three-dimensional fields becomes 

- gV2 (a^^[0,„ e\^] - A,[0^^ e\]])+ ^-g^\<P,,, 4>ki] W\ 0'1+^^7'[03, <P^j\ [03, <P'']+ 

2 

- ^A,7°A* - - /iVI + 2/i03Fi2 . 

(2.18) 

This is the A/" = 8 SYM Lagrangian on IR x 5^ that will be used in computing the ther- 
modynamic partition function of the model. We have cast the contribution of the scalar 
fields (04, . . . , 09) in an SU{'i)ii manifestly covariant form, by rewriting their Lagrangian 
in terms of the 6 representation of SU (4)^, (pij. The spinor fields Aj are four Dirac spinors 
in D = 3 originating from the dimensional reduction of ip. 

Since we will be mainly interested in the finite temperature features of the model, the 
Euclidean version of fl2.18p will be more relevant. It is given by 

C = ^F^.F^-' - 2i\a^D,\' + ^D^0,,D^0*^' + D^<I),D^^<P^+ 

+ gV2 (a^^[0,„ eX'] - A,[0*^ eXj]) + 2tg\[<f)s, A*]+ 

1 1 (2.19) 

- -g'[<P,,, <Pkl] [0^^ 0'1 - -^?'[03, 0i,] [03, <P''] + 
2 

+ ^A,7°A^ + ^0i,-0^^' + /iVs - 2/103^12. 
Z o 

We conclude by noting that, in the above analysis, we have made a particular choice 
in considering the form of the Killing spinor equation. A careful reader might wonder 
if there are other possibilities. Unfortunately, different choices in (12. 4p generally lead to 
inconsistencies: the Killing equation is not integrable or no consistent supersymmetric 
deformation exists. For example, the second type of inconsistency would occur if we 
had simply chosen K^'^ = It is however intriguing to note that the choice K^'^ = 
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(5^ becomes consistent if we alter the background geometry from R x 5^ to AdSs and 
substitute F^^^ with F"^^ or F^^^. In the former case, we would have found a maximally 
super symmetric version of the topologically massive theory, with bosonic symmetry group 
5*0(1,3) X 5*0(7). In the latter we would have instead reached a massive deformation 
of the maximally supersymmetric Yang-Mills with the peculiar interaction Tr((/)3[</)4, ^s]) 
and symmetry group 5*0(1, 3) x 5*0(3) x 5*0(4). This case was already considered in [32] . 
It would be nice to understand better their relations with higher dimensional theories and 
to explore the possible existence of gravitational duals. 



3 BPS vacua and their gravitational duals 

In this section we shall briefly review the structure of the BPS vacua of the A/" = 8 theory 
on R X 5*2 [5] that will be the main ingredients of the thermodynamical investigation of 
section [6] and [71 More specifically, we shall be interested in those vacua that maintain 
both the i?-invariance and the geometrical symmetries. 

In order to have an 5f/(4)/j invariant vacuum, we have to choose (pij = 0. More- 
over, to preserve the invariance under time translations and the 50(3) rotations of the 
background geometry, we require that all the fields are time-independent and that the 
chromo-electric field Ei = F^i vanishes, respectively. The BPS condition can be derived 
from the requirement that on the supersymmetric invariant vacuum the supersymmetry 
variations should vanish. Fermions must be set to zero to saturate the BPS bound and 
consequently the supersymmetry variations of bosons automatically vanish on the vac- 
uum. The supersymmetry variation of fermions, instead, must be set to zero and with 
the above assumptions it reads 

= 5^= [2{Fe^ - - sin ^03)r''^ + 2D^03r^']e , (3.1) 
/i 

[6 and ip are coordinates on 5^) which translates into two simple equations 

Fg^- - sin 9<Ps = 0, D^^3 = 0. (3.2) 
/i 

The reader familiar with YM2 will immediately recognize in these equations, those of 
Yang-Mills theory on the sphere 5*^, for which a complete classification of the solutions 
exists [^B2]. The general solution for a U (N) theory is given by a stack of independent 
U{1) Dirac monopoles of arbitrary charges. In detail, we have 

03 = y Fe^ = -sme A = -r^ [smedip) = -A, (3.3) 

where f is a diagonal matrix with integer entries, for which we shall use the short-hand 
notation 

f={num;n2,N2;...;nk,Nk). (3.4) 
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Each nj represents the Chern-class of the corresponding Dirac monopole and it assumes 
values in Z, while Nj is the number of times that this charge appears on the diagonal. The 
vacuum (13.41) then breaks the original U{N) gauge symmetry to a direct product U{Ni) x 
U{N2) X . . . U{Nk). However, since all fields in (12.181) are in the adjoint representation, 
this breaking will affect the dynamics only through the relative charge (ri/ — nj) between 
different sectors, while the global charge Q = X]/=i Ninj will play no role. 

The gravitational backgrounds dual to the vacua of these theories were derived in [5] 
and further discussed in [61] (where also the relations between vacua of theories with 
SU{2\ 4) symmetry group are studied): they have an S0{3) and an S0{6) symmetry and 
thereby the geometry contains S"^ and factors, the remaining coordinates being time, 
a non-compact variable rj, — oo < rj < oo, and a radial coordinate p. These backgrounds 
are non-singular because the dual theories have a mass gap. The relevant supergravity 
equations can be reduced to a three-dimensional electrostatic problem where p is the 
radius of a charged disk. The ten dimensional metric and the other supergravity fields 
are completely specified in terms of the solution V of the related Laplace equatioio. The 
regularity condition requires that the location where the S"^ shrinks are disks at constant 
T]i (in the p, rj space) while shrinks along the segment of the p = line between two 
nearby disks. The geometry therefore contains three-cycles connecting the shrinking S*^ 
and six-cycles connecting the shrinking S^, supporting respectively non-trivial H'^ and 
fluxes. There is a precise relation between these quantized fluxes and the data of the 
electrostatic problem, namely the electric charges Qi of the disks are related to the RR 
fluxes while the distance (in the rj direction) between two disks bounding a three cycle 
is proportional to the NS flux. To be more specific, this electrostatic description of a 
non-trivial vacuum generically contains k disks, whose positions are parameterized by k 
integers n/ through the relations 

m = (3.5) 

These integers are identified with the monopole charges nj in (13.41) . Moreover each disk 
carries a charge Qj given by 

Qi = — (3.6) 

where Nj are the same integer numbers counting the degeneracy of each monopole charge 
in the gauge theory. At the level of supergravity data, the above picture realizes k groups 
of D2 branes, each of Nj elements, wrapping different two-spheres. This is the geometric 
manifestation of the breaking of the gauge symmetry to a direct product U{Ni) x U {N2) x 
■ ■ ■ X U{Nk). The charges nj instead combine into NS5-fluxes given by rij — nj. Again 
the total charge seems to play no role. 

In our field theoretical analysis we have neglected the time component of the gauge 
field Aq, which disappears from (13.21) when considering the solutions (13. 3p . Its dynamics 
is implicitly governed by the requirement that Ei = 0, which, for a time-independent 
background, becomes DiAo = 0. It is a trivial exercise to show that the most general 



•^This problem has been recently tackled in [13] and [33], searching for a dual description of Little 
String theory on 
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solution of this equation is provided by = when the topology of the time direction 
is R. In the finite temperature case where time is compactified to a circle 5*^, the most 
general solution is, instead, given by Aq = a, where a is a constant diagonal matrix, 
namely a flat-connection living on 5*^. This will play a fundamental role in studying the 
thermodynamical properties of the theory. 

It is instructive to look at the BPS vacua also at the level of the Euclidean Lagrangian: 
this will elucidate the emerging of an interesting class of instanton solutions thoroughly 
studied in [15]. If we focus on the bosonic sector of our model and we set (pij = to 
preserve the SU{4)j^ symmetry, we can write 

= ^F^pF'^P + ^D^<P3D^(f)3 + V9f^'<Pl - 2fihFg^. (3.7) 

This Lagrangian can be easily arranged in a BPS-form, i.e. as a sum of squares and total 
divergences. In fact, after some algebraic manipulation, the Euclidean Lagrangian can be 
cast in the following form 

= ±- smeDt{<Pl) T /^„(03i^/3pe"^O + sin^ ( F^e ± -^D^iP^] + 
^ V sm^ ^ y 

1 2/1 
+ {Ft^ T sin dDe(l):,f + [Fe^-— sin ^(/i03 t Dt^3 
smt/ smO \ jj,"^ 

Consequently, the minimum of the action is reached when the fields satisfy the following 
BPS-equations 

(a) : Fte±^^ = (6) : Ft^TsinODeh = (c) : Fe^-\ sin 6 {fi(f),TDth) = 0, (3.9) 
smt^ yLi^ 

or in a compact and covariant notation 

VgepuxF''^ = t2/^p03 + 2/iA;^03, (3.10) 

where kp is the Euclidean version of the time-like Killing vector of the metric on R x 5*2. 
The vacuum equations (13. 2p are just a particular case of (13. 9p or equivalently (13.101) . They 
emerge when we add the requirement of time-independence and vanishing of the chromo- 
electric field Ei. From (13. 8p it is manifest that all our vacua (13. 3p possess a vanishing 
action and they are all equivalent from an energetic point of view. 

It is natural to ask now what is the meaning of the Euclidean time- dependent solutions 
of ( 13. 9p . The action on these solutions reduces to 

Sciass = T- [ dOdip sin e [ dtdtTT{(j)l) , (3.11) 

which is finite, and thus relevant for a semiclassical analysis of the theory, if and only 
if 03(t = — oo) = and 03(t = oo) = ^f^- In other words, these solutions are 
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interesting if and only if they interpolate between two vacua: one at t = — oo and the 
other at t = +00. Their finite action is then given by 

Sciass = T- / sm9d9d^ / dtd,TT{4>l) = T^-^mfJ " MtJ) , (3-12) 

where we have reintroduced the relevant coupling constant factors. We recognize the char- 
acteristics of instantons in these (Euclidean) time-dependent solutions. At the quantum 
level, they will possibly induce a tunneling process between the different vacua. At zero 
temperature Lin [15] discussed the effect of these instantons from the gauge theoretical 
side, at weak coupling, and from the gravity side, that should describe the strong-coupling 
limit of the theory (see also [IS]), finding precise agreement in both regimes. Moreover 
he argued, in analogy with the plane-wave matrix model, that because of the presence 
of fermionic zero-mode^ around these instanton solutions, the path-integral for the tun- 
neling amplitude is zero. The vacuum energies would not be corrected and the vacua 
are exactly protected at the quantum mechanical level: in particular they should remain 
degenerate. This kind of instantons has also been recently considered in [47j . 

In the rest of the paper, in any case, we shall neglect the effect of these solutions 
since we shall work at zero-coupling and in this limit the probability of tunneling is 
exponentially suppressed anyway. 



4 Free SYM partition functions in monopole vacua 

In this section we shall derive the finite temperature partition function in the BPS vacua 
(13. 3p . taking the limit R —>■ 0. We follow a path-integral approach where the compu- 
tation is reduced to the evaluation of one-loop functional determinants in the monopole 
backgrounds. Since at finite temperature the Euclidean time is a circle 5*^ of length 
f3 = 1/T, we can also allow for a flat-connection a wrapping this S^. The mode a will 
play a very special role because it is the only zero-mode in the decomposition into Kaluza- 
Klein modes on S"^ x S^. Consequently, as stressed in [9], the fluctuations described by a 
are always strongly coupled, including in the limit i? — > 0. 

When the vacuum is trivial, there is no breaking of the U{N) gauge symmetry and 
the final result for the partition function is given by a matrix integral over the unitary 
matrix U = exp [z/3a] 

Z{f3) = j [dU] exp 1^ ^ [^b(x") + {-IY+^ zpix"")] Tr(f/'^)Tr(f/-")| . (4.1) 

The functions zb,f{x) are respectively the bosonic and fermionic single-particle partition 
functions (here x = e~^), counting the one-particle states of the theory without the 

^The instantons are 1/2 BPS solutions and therefore we expect 8 fermionic zero- modes associated to 
the broken supersymmetries 
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degeneracy coming from the dimension of the representation (the adjoint representation 
Adj in our case) and without any gauge invariant constraint 

.^,Hx) = $:e-/^^l""\ (4.2) 

i 

The exphcit form of the thermal partition function is obtained by integrating over the 
matrix t/ [H [9] 

oo oo oo oo 

ni=0 n2=0 mi=0 m2=0 

# of singlets in {sym"^ (^c(7) ® sym"^ (Adj) ® • ■ • 

(g) antisym""^ (Arfj) (g) antisym'"^ (Adj) ® ■ ■ ■ } : (4.3) 

the partition function is expressed as a sum over the occupation numbers of all modes, 
with a Boltzmann factor corresponding to the total energy, and a numerical factor that 
counts the number of singlets in the corresponding product of representations. Particle 
statistics requires to symmetrize (antysimmetrize) the representations corresponding to 
identical bosonic (fermionic) modes. 

The same result can also be obtained starting from 

Z{(3) = Tr [e-^^] = Tr [x^] , (4.4) 

where H is the Hamiltonian of the theory. To calculate (14. 4p at zero coupling we need a 
complete basis of states of the free theory or, thanks to the state-operator correspondence, 
of gauge-invariant operators and we should count them weighted by x to the power of 
their energy. A complete basis for arbitrary gauge-invariant operators follows naturally 
after we specify a complete basis of single-trace operators. At the end, one can write (14.41) 
in terms of single-particle partition functions z§ p{x) [9] as 

Z{P) = [ [dU] exp I 5^ f; - [z^ixn + i-ir^'z^ix-)] XRiUn ] , (4.5) 

[ R n=l^ J 

where the sum is taken over the representations R of the U{N) gauge group@ and Xr{U) 
is the character for the representation R. The result (14. ip is reproduced when all fields are 
in the adjoint representation: the variable U has to be identified as the holonomy matrix 
along the thermal circle, i.e. the Polyakov loop. The path-integral approach provides 
therefore a physical interpretation for the unitary matrix U, otherwise missing in the 
Hamiltonian formalism. On the other hand the Hamiltonian construction explains how 
the group integration forces the projection into color singlets and how it emerges the 
structure of the full Hilbert space. 

^We consider the possibility to have fields in an arbitrary representation. 
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From the previous results we learn that once the representation content is specified, 
the full partition function is completely encoded into the single-particle partition func- 
tions z§ p. However, the structure of the gauge group is more complicated on monopole 
backgrounds, consisting into a direct product of U{Nj) factors: consequently our con- 
stituents fields transform also under bifundamental representations, producing additional 
complications for the explicit expression of the matrix model. We also remark that bifun- 
damental fields can transform non-trivially under U{1) rotations and implementing the 
Gauss's law hides some subtleties in three dimensions, when background monopole fluxes 
are present [23] : this potential additional freedom could affect non-trivially the spectrum 
of physical operators in our theory. For the theory we are investigating, however, the 
free-field spectrum is simply obtained by truncating the four-dimensional parent theory, 
suggesting that the A/" = 8 counting is conveniently performed through the relevant U{1) 
projection on the A/" = 4 single-particle partition functions. This is what we do in ap- 
pendix C, where we construct the projector that eliminates all the fields which are not 
invariant under the U (1) and we derive, even in the non-trivial vacuum, the single-particle 
partition functions for bosons and fermions. While this is certainly the quickest way to 
obtain these quantities, we prefer to adopt here a path integral approach which in turn 
provides also the contributions of fermions and bosons to the Casimir energy and allows 
for a careful treatment of the fermion zero modes. In the path-integral computation all 
the subtleties will be treated in the well-defined framework of the (^-function regularization 
procedure and in this section we present only the final results, referring for the technical 
details to appendix B. 

4.1 Scalars 

Let us first describe the contribution of the six SU (4)/j scalars (pij to the partition function 
in the background (13.31) and in presence of the fiat-connection a: it amounts to the eval- 
uation of the determinant of the scalar kinetic operator. We have to solve the associated 
eigenvalue problem, i.e. 



where the hatted quantities are computed in the relevant background. In the following 
we shall drop the subscript ij and we shall consider just one field denoted by 0. The 
total result at the level of free energy is then obtained by multiplying by six the single- 
component contributions. Since is a matrix-valued field, we shall expand it in the 
Weyl-basis, whose elements are the generators Hi of the Cartan subalgebra and the ladder 
operators E" 




(4.6) 




(4.7) 
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We shall also expand the background fields in this basis and define the following two 
accessory quantities 

= {a\a) and qa = ^ ■ (4.8) 

Here denotes the projection of the fiat-connection a along the root a and is the 
effective monopole charge measured along the same root. Once the time-dependence 
is factored out, the original eigenvalue problem splits into two subfamilies: A^(A^ — 1) 
independent eigenvalues coming from each direction along the ladder generator and — 1 
independent eigenvalues coming from the directions along the Cartan subalgebra. We 
can simply focus our attention on the first family, since the latter can be obtained as a 
limiting case for Cq,, g^, — > 0. The relevant eigenvalue equation can be solved algebraically 
if we introduce the angular momentum operator in the presence of a U{1) monopole of 
charge g^, as explained in appendix B, and the resulting spectrum does not depend on 
the sign of g^. By using (^-function regularization, the scalar contribution to the effective 
action can be easily computed as 

(\n I °° rascal. ( n\ \ oo scal.(^n\ 

^ m - 1) + E ^^-^e-^- +(iV- 1) (4.9) 
n=l / ri=l 

where the scalar single-particle partition function is given by 

.,7^-(x)=xl-l~^V2^_I±|_ + ^y (4.10) 

4.2 Vectors 

Evaluating the contribution of the system (A^, ^3) is more subtle and involved: the fields 
are coupled through the Chern-Simons term and the Lagrangian for requires a gauge- 
fixing procedure, with the consequent addition of a ghost sector. A convenient choice for 
such a gauge-fixing appears to be 

£3./, = (^,^^^-^[03,03])', (4.11) 

where 03 = y and the hatted derivative is defined in fIB.lQp . With this choice some 
of the mixing-terms in the Euclidean quadratic Lagrangian cancel and we obtain the 
relevant eigenvalue-problem for computing the vector-scalar contribution to the partition 
function: it is defined by the system of coupled equations, written explicitly in (]B.26p . 
Since both the geometrical and the gauge background are static, the time-component 
of the vector field Aq decouples completely from the eigenvalue system and satisfies the 
massless version of the scalar equation previously studied. For the moment we shall forget 
about Aq since its contribution will be cancelled by the ghost determinant. We are left with 
a purely two-dimensional system where all the indices run only over space: the spectrum 
is again conveniently determined by factoring out the time-dependence and projecting the 
eigenvalue equations on the Weyl basis. We remark that the equations involve also the 
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Laplacian on vectors in the background of a monopole of charge g^, besides the Laplacian 
on scalars. The full computation of the spectrum is reported in appendix B: we obtained 
three families of eigenvalues, denoted by A+, A_ and A3. The contribution of A3 will be 
cancelled by the ghost determinant and we just consider, at the moment, the first two 
families A±, which instead yield the actual vector determinant in the roots sector 



{K + 5g.) - 2 X: ^^4^^ e-^- , (4.12) 

n=l / 



where 



zZ-{x)=x 



Ax , ^ 1 + a; 
2q, 



(4.13) 



1 — xy 1 — X 

We remark that the results ( 14.12^ and ( 14.130 were shown to hold under the initial as- 
sumption > 1. The extra-cases to be considered are = 0. By recomputing the 
spectrum for = 1/2 we get the same results: quite surprisingly this does not happen, 
instead, for Qa = and we get 

rr(ga = 0) = -2V^^^^e^"^'^'^ with z-=(x) = -^, (4.14) 

n=l ^ ' 

a factor —1 missing in the limit. To complete the discussion, we notice that, when 
multiplied by (A^ ~ 1)? (14.141) is the contribution of the Cartan components; the results 
(14.121) and (I4.13P extends also to negative charges by simply replacing g^ with Ig^]. 



4.3 Ghosts and ^0 

Let us discuss now the contributions to the partition function of the eigenvalues A3, of 
the field and of the determinant of ghost operator 

-□■+[03, [03,-]] : (4.15) 

they do not cancel completely but, importantly, they give a measure of integration for the 
flat- connect ion. It is possible to show that when g^ 7^ we have a complete cancellation 
of the different contributions: crucially for g^ = this does not happen and a modification 
of the measure for the fiat-connection is induced 

n 2,e-%'rin(^)= n ''-"'(^)- (''■16) 

a G roots ^ ^ a £ positive roots ^ ' 

with ga— with ga — 

The meaning of this measure is quite transparent: the monopole background breaks the 
original U{N) invariance to the subgroup 117=1 ^(^/)' (I4-16P being the product of the 
Haar measure of each U {Nj) component, as can be easily checked by recalling the explicit 
form of the roots and the definition of go,. As a matter of fact, in non-trivial monopole 
backgrounds, when we shall write the integral over the fiat-connections we will be naturally 
led to consider a unitary multi-matrix model instead of an ordinary one. 
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4.4 Fermions 



The contribution of the fermions to the total partition function needs a careful analysis. 
At first sight, apart from having antiperiodic boundary conditions along the time circle, 
the computation of the fermion determinants seems to follow closely the bosonic cases. 
We have again A^(A^ — 1) independent eigenvalues coming from each direction along the 
ladder generators and TV — 1 independent eigenvalues coming from the directions along the 
Cartan subalgebra, that can obtained as limit of vanishing flux. The computation of the 
spectrum is quite technical as in the vector case and boils down in solving the eigenvalue 
problem for a family of effective massless Dirac operators 2)*^"^ (see app. IB.4p on the two- 
sphere, in the effective monopole backgrounds provided by g^. The spectrum of D^"\ as 
expected in two dimensions, consists in a set non-vanishing eigenvalues, symmetric with 
respect the zero, and in a finite kernel, as predicted by the Atiyah-Singer theorem. These 
zero-modes are chiral and can be classified by using the eigenvalues of the operator (cx-f), 
playing the role of 75: we shall denote z/± the number of zero modes with eigenvalue ±1. 
A simple application of the index theorem shows that = \qa\—(la and i/_ = \qa\ + qa, 
namely for positive we have only zero modes with negative chirality and viceversa. As 
shown in appendix IB.4t the contribution of the first set of eigenvalues to the effective 
action can be easily evaluated 

rf =E f-f (2|^"l'+ 3k.P+ M-f: ^<!r-(^")e^^""" 1, (4.17) 

aGrootsV n=l 



with 



1 1 , 



Next we consider the contribution of the zero-modes of the effective Dirac operators: 
in a monopole background, this subsector originates the spectral asymmetry [18] of the 
three dimensional fermionic operator and therefore the potential appearance of a parity 
violating part in the effective action. In particular, we could expect the generation of the 
Chern-Simons anomalous term (we refer to [331 [M] fo'^ a complete discussion of this issue). 
Concretely, in our case, the explicit computation of the zero-mode contribution amounts 
to evaluate a family of one-dimensional massive fermion determinants, in a flat-connection 
background (see appendix lB.4p . It is well-known that the (^-function regularization scheme 



carries an intrinsic regularization ambiguitj{§ in this case, depending on the choice of some 
branch-cuts in the s-plane, affecting the local terms in the effective action [531 El]- For 
us all the different possibilities boil down to two alternatives: we can regularize the 
contributions associated to the zero-modes of negative and positive chirality by choosing 
opposite cuts in defining the complex power of the eigenvalues (one on the real positive 
axis and the other on the real negative axis) or by choosing the same cut. We find quite 



^This ambiguity is not something peculiar of the C-function regularization, but it appears in different 
forms also in other regularizations: in the usual Pauli-Villars approach, for example, this ambiguity 
translates into a dependence of the local terms in the effective action on the sign of the mass of the 
regulator. 
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natural to use the same procedure for all the four fermions present in the theory: we surely 
preserve the i?-symmetry and the global non-abelian symmetry in this way. Within this 
choice, the following results hold from our one-dimensional fermion determinants: taking 
opposite cuts we get 

^Ia = E (1 - ^)f^^^ (^^ + ^) - E E ^2|g.|x"l-le^^-xt. (4.19) 

«groots ^ ' agroots n=l 

Here r = ±1 and its specific value depends on the cut selected for the zero-modes of 
positive chirality. Choosing instead the same cuts we obtain 

^Ib= E [/5/i(k«r + ^)+^r/3a„g„]- f2t^2\q^\x-\'^-\e'^-^-x^^. (4.20) 

aSroots ^ ' aSroots n=l 

Again r = ±1 according to the specific choice of the cut (real positive or negative axis): 
we must stress, however, that this ambiguity will become irrelevant when we shall perform 
the integration over the flat-connections. 

We remark that there is an important difference between the two expressions: in the 
second case we have a new term in the effective action, depending explicitly on the fiat 
connection. To understand its nature, it can be equivalently written as 



^/3g„a„ = irp{NTT{af) - Tr(a)Tr(f)). (4.21) 



We immediately recognize the SU{N) part of the usual Chern-Simons term, calculated 
in our particular background. The related regularization choice is therefore consistent 
with the intrinsic parity anomaly of three dimensional gauge theories. We stress that the 
above contribution arises just in the monopole vacua and it is related to non-perturbative 
properties of the fermion determinants. We also observe that the two results differ in the 
charge- dependent contribution linear in /?, and we will see this to modify crucially the 
Casimir energy. 

Summing now, in both cases, the kernel contribution to Ff we get 

ri =E (-ff 12r|g.p+ (3r + l)\q^\) -f^ t^z^/J-ix-)e^^-A, (4.22) 

«groots\ n=l / 



with the first choice and 




oo 



+ irq^a^ ^^^^r' (x")e^^"'^'^ ), (4.23) 



n 

n=l 

in the latter. Happily the single-particle partition function is the same for both the 
regularization choices 

zZ-ix) = 0:1-1 + (xi + x-^). (4.24) 

The contribution of the Cartan components is of course obtained from the above results 
by simply setting = 0. 
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4.5 Partition functions 



The next step is to collect the different contributions, coming from the functional deter- 
minants, and write down the total result as a compact integral over unitary matrices. 
According to the previous discussion, we must distinguish two cases, depending on the 
form of the spinor determinant fl4.22p or fl4.23p . We shall first consider the choice fl4.22p . 
The complete effective action, obtained by including roots and Cartan contributions with 
the appropriate multiplicities, can be expressed as 



Seff. =-PVo+ E + <r(^") + (-l)"+^4^,^r-(x"))e^'^^- + 

oSroots n=l 
oo ^ 

+ (iV - 1) y -{Gzr'ix'') + z'o'-i^l + (-l)"+'42r"-(x")) = (4.25) 

n=l 

oo oo 

^-PVo+ (x")e*"^'^'» + (iV-l)^-4°*-(x"), 



n " ^ — ' n 

aGroots n=l n=l 



where we have introduced the total single-particle partition functions and the Casimir 
energy Vq of the configuration 

Vo = r J2 (Mlal' + M- (4-26) 

agroots 

The matrix structure hidden in (I4.25P appears manifest when writing the original Polyakov 
loop U = exp(i/3a), associated to the diagonal flat-connection a, through k sub- matrices 
Uj acting on the invariant subspaces implicitly defined by the monopole background (13. 4p . 

The Nj X Nj unitary matrices Ui have the form Uj = diag(e*'^"i, . . . ,e^°''^'), where we 
have parameterized the original fiat connection a as follows: 

a = diag(a}, . . ■ ,a]v^, Qi, ■ • • , a^,^ , , a(, . . . ,a^^, • • •)• (4.27) 

Let us consider now the subset Ajj of the positive root^ of SU (N) whose first and second 
non vanishing entries belong respectively to the and J*'* invariant subspace of f. The 
effective charges Qa = = and, consequently, the z*°^' take always the same value 

for this class of roots. The sum over roots on this subset reduces to 

Ni Nj 

J2 e'"^"" = e*"^^"^-"/^ = Tr([/7)Tr(?7j"); (4.28) 

aeAij i=l j=l 

the analogous subsector Ajj given by the negative roots yields Tr(f/j")Tr(f/j). We remark 
that the pre-factor zi°*" is however the same for both cases since it depends just on the 



^ The roots of SU{N) are all the N{N - 1) permutations of the iV-vector (1, -1, 0, • • • , 0) and they 
can be separated in positive and negative according to the sign of the first non zero entry. 
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modulus of the effective monopole charge. The subset of roots Bj whose first and second 
non vanishing entries hve in the same J*^ invariant subspace of f have instead effective 
monopole charge zero. Then the contribution of this subsector is simply given by 



n=l aSi?/ n=l *7^i=l 



Ni 

inf3{al-a'.) 



oo 



E 



° ^ ^(Tr(f/f)Tr(f/7) - AT^) 



(4.29) 



n 

n=l 



Because of the results f l4.28p and (14.291) . it is convenient to change our notation and 
to define the k x k matrix- valued single-particle partition function zf'j': the diagonal 
elements •, the off-diagonal ones are instead identified with the function 

associated to the charge The matrix zf'j' is symmetric since everything depends 

just on the modulus of the charge. The complete effective acton takes the elegant form 

oo oo 

Seff. = -PVo + (^")Tr(f/7)Tr(f/J") - -zfH^l- (4-30) 

IJ 71=1 71=1 

The last term drops if we consider U{N) instead of SU{N). Remarkably the structure 
of the matrix action is perfectly consistent with the measure found in (14.161) . which is 
exactly the Haar measure for this multi-matrix model. 

The above analysis is practically unaltered when considering the fermionic contribution 
(14.231) in the effective action, except on a couple of points. It changes the value of the 
Casimir energy Vq, which now vanishes identically, and we have a new important addition 
to (I4.30p . that can expressed in terms of the determinants of the unitary matrices Uj 

ndet(?7,)^(^"^-«M =r5^(iVn,-g)log(det(f/,)), (4.31) 

7=1 / 1=1 

where Q = X]/=i NjUi. As a first remark, we notice that new contributions depends still 
on the differences rij — nj, consistently with the decoupling of the total U{1) charge of the 
monopole configuration. Then we observe that the two different values r = ±1, related 
to our regularization choice, produce the same result when integrating over the unitary 
group: the difference can be reabsorbed just changing integration variable Uj ^ (f//)"^, 
which leaves the measure and (14.301) unaltered. From now on, we shall set r = 1. 

In the trivial vacuum we obtain a partition function that is a straightforward gener- 
alization of the unitary matrix model discussed in [9] 

Z = y"rff/exp f ^^z*°*-(x")Tr(f/'^)Tr(f/t") j (4.32) 



.n=l 



where the function z^^'{x'^) encodes the dynamical content of the three-dimensional su- 
persymmetric theory. Notice that the Casimir energy is identically zero, since it vanishes 
for each contribution both bosonic and fermionic. 
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The situation changes in non-trivial monopole vacua: we get respectively 

2a = [ f[[dUj] exp (-pVo + Y.fl (^")Tr(f/7)Tr(f/J") ] (4.33) 

1=1 \ IJ n=l J 

and 

/k / oo ^ \ ''^ 

Y[[dUi] exp ^ -zfj-(a;")Tr(f/7)Tr(f/J") JJ det(t/j)(^"^-«), (4.34) 

1=1 \ I J n=l ^ / 1=1 

depending on our regularization choice. First of all we see that the partition function 
is related to a unitary multi-matrix model: the gauge group is broken in factors and 
states in the bifundamental representation are present, with energies clearly encoded into 
the off-diagonal entries of the single-particle partition function z^'j' . Let us discuss on 
general grounds the effects of the different choices for the fermion determinants. A first 
mild diversity arises in the Casimir energies: from (I4.22p we have a non- vanishing Vq, 
with arbitrary sign, while (14.231) leads to a vanishing result. We recall that the Casimir 
energy is supposed to correspond to the mass of the dual geometry [39]: in the first 
case it seems that different backgrounds supports different, monopole dependent, masses, 
suggesting a possible lifting of the vacua degeneracy at quantum level. The second choice 
is instead consistent with the believed degeneracy: unfortunately no computation from 
the gravitational side seems to be available up to now and we do not have further insights 
on the meaning of the different results. 

The presence of the new terms (I4.3ip in the matrix model (14.341) can be, instead, better 
understood at the level of partition functions. First of all we notice that the matrix 
integral implementing the Gauss's law is actually over unitary matrices t//: the U{1) 
phases contained into the the f//'s play a non-trivial role in the monopole background. 
This has to be contrasted with the trivial vacuum: there the effective action is invariant 
under U{1) transformations and we can simply forget the integration over the center. In 
the non-trivial vacuum the resulting effective action (I4.30p is not invariant under phase 
rotations, as an effect of the off-diagonal terms in the single-particle partition function, 
and the U{1) integrations precisely correspond to selection rules in the bifundamental 
sector. It is not difficult to realize that within the first regularization the matrix integrals 
select states having vanishing U{1) charge, with respect to all U{Nj) group factors. To 
understand the effect of the new terms in (I4.34p instead, we simply observe that the 
determinants depend just on the U{1) phases and modify non-trivially the selection rules 
of the bifundamental sectors, according to the charges of the monopole background. We 
shall say in this case that our regularization procedure correspond to the choice of a 
charged vacuum, as discussed in [26], while we will refer to the first possibility as to 
the uncharged vacuum. Since at the quantum field theory level both choices seems to 
be allowed, we think it is instructive to investigate the thermodynamics in both cases, 
deferring a deeper understanding of the different possibilities to future studies, in the 
context of supersymmetry and gravitational duals. 
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We end this section introducing the simple modification to the effective action due to 
chemical potentials for the SU{4) i?-charge. In the path integral approach their effect 
amounts to simply adding an imaginary S'f/(4) flat connection = i{QiQY- + Q2Q^ + 
Q^Qf') in the Euclidean time direction. Here Qf^ are the Cartan generators of SU{A) and 
R denotes the relevant representation: 4 for the spinors and 6 for the scalars. One flnds 
the new partition functions 

/ A 

2x 

X) 



V ^ ^ p=i 

p=i / 

with y = and 

=^{-ni + n2-ns) =^{-^1-^2 + ^;). (4.36) 
5 Thermodynamics in the trivial vacuum 

We have seen in the previous section that the thermodynamics in the trivial vacuum 
is governed, in the zero-coupling approximation, by the one-component unitary matrix 
model 

Z = y^t/exp ^^^4''*-(x")Tr(f/")Tr(f/^")j (5.1) 

where the function Zq°^'{x'^) encodes the dynamical content of the three-dimensional su- 
persymmetric theory. Notice that the Casimir energy is identically zero, since it vanishes 
for each contribution both bosonic and fermionic. 

When is large we can trade the integration in (15.11) over the unitary group for an 
integration over the normalized distribution function p{6) of the continuous eigenvalues 
e*^ of U, with —IT < 6 < 71. More precisely we can write the integral over the unitary 
matrices in terms of the Fourier-modes (p„,p„) deflned as 

^ oo 

p{0) = - + J2(Pn^''''+Pne-'''') ■ (5-2) 
n=l 

Following [H [9], we can then reduce the integral to the standard form 

with V^(x") = -(1 - 4°*-(x")) . (5.3) 
n 
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In the large hmit, (15.31) is dominated by the absolute minimum of the quadratic action 
S = ^ PnPnV{x^) which is reached for Pn = for every n if y(x'^) is positive definite. For 
small temperatures, namely small x, the function V{x^) is positive for any n and close 
to 1/n since l^(x") ~ i for X < 1. (Recall that Zq°^-{x'^) vanishes as x approaches zero.) 
Therefore the partition function is 1 at the leading order and it is simply given by the 
small fluctuation around the minimum at the subleading order: 

ZocITt (5.4) 

n=l ^ u \ // 

When we increase the temperature, x approaches 1 and the above description is reliable 
up to the smallest value Xc where V{x"') becomes negative. Since zl"*'{x) is a monotonic 
function ranging from to infinity, this value always exists and it is reached for n = 1, 
namely 

V{x,) = l-4"\x,) = 0. (5.5) 
This algebraic condition, whose explicit form is 

Vix,) = 1 - 4°'- = 1 - iAzf""- + 6^o""' + 4"-) = 

_ + l)^ {xc - ^Xc^l^ + - 4^ + l) _ ^^-^^ 

^ (T^^^J^ 

can be exactly solved, since it can be reduced to an equation of fourth degree. It possesses 
just one solution in the interval [0, 1] given by 

Xe = (^2 + 2V2 - + 8^2^ ~ (0.104688)1 (5.7) 

It is interesting to compare this value with the critical temperature computed in [39] for 
A/" = 4 on ^'^/Zfc. This theory should in fact reproduce our model when k goes to infinity. 
However, the three-dimensional theory obtained in this limit lives on a S'^ sphere whose 
radius is half of the radius of the original S"^: this means that Xc = limk^oo xl{k) . To 
facilitate the comparison with the four dimensional literature and in particular with the 
results of [201 EH] in what follows we shall replace the basic variable x with x^. In [39] the 
Xc{k) for = 10 is 0.104689 which is already very close to (15. 7p . 

Above this critical value the integral (15.31) is no longer dominated by the trivial min- 
imum Pn = Pn = and one has to look for other saddle-points [H [9]. Following [9], one 
can easily show that above Xc the dynamics is governed by a distribution different from 
zero only in the interval [—6*0,6*0] and given, in first approximation!^, by 




PW = -3TiYJ-nM::^ -sinM^ with cosM^ UJl--^. (5.^ 



7rsin^(f)V"" V2y \2J ^ 4°'ix) 

This behavior at Xc produces a first-order transition with the same qualitative character- 
istics of the four-dimensional model. 



^We are assuming that the relevant features are completely captured by the first mode n — 1. 
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5.1 Chemical potentials 

A natural and intriguing generalization is to add chemical potentials for the /^-charges, 
while maintaining the trivial vacuum as a gauge background. 

The critical equation has still the form (15.61) but iz^^^^' and 6zQ'^°^'^^ are substituted by 



4^, 



6z, 



spin. 



2x 



'1-x) 



p=i 



ical. 



^ ' p=l 



(5.9) 



which is (I4.35P for = 0. The effect of small chemical potentials can be easily computed 
by treating them as a perturbation and expanding around (r2i,r22;^3) = (0,0,0). This 
yields the following result 



Tnin) = Th{0) - 0.113946 ^ - 0.054438 fj Qi 

1=1 i=l 

3 

- 0.036442 0.014059 ^ QfQ'^j + 0{Q^ 



(5.10) 



1=1 



i<j 



where all the numerical coefficients are actually known exactly, but their explicit ex- 
pression is long and irrelevant. The presence of small chemical potentials decreases the 
Hagedorn temperature. 

In figH] we display the dependence of the critical temperature Th for the three par- 
ticular choices of critical potential (f2i,f22,^3) = (f^,0, 0), (f2i,f22,^3) = (f^, i^, 0) and 
(fii, f22, ^^3) = (f2,f2,r2). In all three cases, the behavior around f2 = 1, in a trivial vac- 
uum background, is similar to that of the J\f = A theory in four dimensions discussed 
in [I9l[20]. We find, in fact: 



(l]i,fi2,f^3) = (^^,0,0) : Th = 



1 



iog(i - ^) L 

l-Vt 



log 2 

l-Vt 
log 4 



1 - 



• log(-log[l-l]]) 

iog(i -n) 

g 2(i-n) _|_ (9(g 



log 2 
30 
log 4 



(5.11) 



log 4 



e ^(i-n) + (9(e •^(i-n)) 



which have the same qualitative behavior of the analogous equations found in [20] for 
the A/" = 4 theory. This similarity suggests the possibility to consider decoupling limits 
analogous to those performed in [20] for the J\f = 4 theory. This might help to single 
out some subsectors of the present model with simple properties at the (full) quantum 
level [221 [231 [21] • However, this analysis is left for future research. 
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Figure 1: The continuous, dashed and dot-dashed hues correspond to (fii, r22, i^a) = (51,0,0), 
(ill, £72, ila) = (fljfl, 0) and (fii, i72, f^s) = (ri,fi, il) respectively. Ah the curves reach = 1 when 
Xc approaches zero. 



5.2 High temperatures 

In the high temperature regime the eigenvalue distribution becomes almost like a delta- 
function [9]. Therefore p„ = 1 and the free energy can be evaluated by looking at the 
expression of the functional determinants in the background of vanishing flat-connections. 
When the chemical potentials are strictly zero the leading contribution to the free energy 
F = -ThgZ is (see (ES]), flKl4D and (IRTSD ) 

F = --Ci3)V{S^)N^T'^ + 0{T^). (5.12) 

71 

We see that the limiting free energy density here coincides precisely with that of the 
A/" = 8 super Yang-Mills theory in fiat three-dimensional space. Taking the dimension- 
less parameter TR to infinity is equivalent to taking the limit of large volume at fixed 
temperature, loosing in this way any memory of the original deformed supersymmetry. 
We can also notice that no dependence appears, at the leading order, on the particular 
monopole vacuum on which the expansion has been performed and the result fl5.12p is 
actually general. 

It is interesting to consider the corrections to this result when chemical potentials 
are taken into account. The first non-trivial contribution is easily evaluated by using the 
expansions of Lis (2;) presented in flB.lSp : we simply notice that chemical potentials appear 
as imaginary parts of the flat-connections and are contained in the variable z introduced in 
the appendix lB.il Summing carefully the contributions coming from bosons and fermions, 
we obtain the free energy 



F = -V{S'')N^T^ 



Ic(3) + t|(3->o4 



0{T^), (5.13) 
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where we introduced the relevant combination Ui = Qi/T. This resuh is perfectly con- 
sistent with the computation performed in [l9], for a system of free D2 branes in the 
presence of chemical potentials. 



6 Thermodynamics in non-trivial vacua I 

We shall first consider the multi-matrix model, (14.331) . which originates from the uncharged 
vacuum. We recall that in this case the partition function is defined by the matrix integral 



2, 

where 



k / CXD \ 

l[[dUj] exp -f3Vo + (^")Tr(t/7)Tr(t/t") , (6.1) 

1=1 \ IJ n=l J 



Vo = rJ2 (4ka|' + |g.|) (6.2) 

agroots 

is the Casimir energy. The value of the Casimir energy is puzzling not only for the r 
dependence, making its sign ambiguous, but also because it depends on the charge of the 
vacuum so that it is different for different vacua. At the supergravity level we expect 
instead these vacua to be degenerate. This last feature is reproduced within our second 
regularization choice, giving a vanishing Casimir energy and consequently degenerate 
vacua: the price we pay is the introduction of the logarithmic interactions (I4.34p that will 
be studied in the next section. 

The large A^-limit of the matrix-model (16. ip is investigated by generalizing to a multi- 
dimensional case the technique presented in the previous section: we introduce the density 
functions pi{Oj) associated to the matrices Uj and in terms of the Fourier- modes p/„, 

Pi{Oi) = — + V(p/ne^"^^ + P7„e-^"^0, (6.3) 
27r ^ — i 

n=l 

the matrix integral (16. 3p reduces as well to an infinite set of independent gaussian integrals 



„ A, ^ ^ 

2a= Y\. DpinDpIn exp -/3 Vo - ^ ^ PlnpJn " (5/ J 
7=1 IJ n=l ^ 



(6.4) 



ij 11— ± V 

\ Vij{x") I 



where we have introduced the filling fractions s/ = Ni/N. In the large A^ limit (16.41) is 
dominated by the absolute minimum of the quadratic action 

oo 

5 = 5^5^P/nPJnV^7j(x"), (6.5) 
IJ n=l 

which is given by p/„ = for every I and n if the quadratic form Vij{x^) is positive definite. 
For small temperatures, namely small x, the eigenvalues of the matrix Vij{x"') are all 



28 



positive and close to 1/n since Vui^x"^) ~ ^6ij for x ^ 1 (we recall that zf'j'^x"') vanishes 
as X approaches zero). Therefore the partition function is simply given by the Casimir 
contribution at the leading order and by the small fluctuation around the minimum at 
the subleading order 

Za oc e-^^° n ^ , (6.6) 

AAdet(V,j(x-)) ^ ^ 

When we increase the temperature, x approaches 1 and the above description is reliable 
until the quadratic form V/-j(x'^) develops the first negative eigenvalue. This occurs at the 
smallest Xc for which one of the eigenvalues of V/j(x") vanishes, or equivalently for which 

det {Vjj{x:)) = 0. (6.7) 

The smallest Xc = e"^/^'", namely the smallest critical temperature, is obviously obtained 
for n = 1 which provides the strongest condition. Moreover this critical value always 
exists since zf^j'{x) is a monotonic function ranging from to infinity when x G [0, 1]. 

We are now ready to investigate the dependence of the critical temperature on the 
non trivial monopole background. We start by considering a configuration f with just two 
sectors of equal length. It is given by 

f = (ni, ... ,77,1,77,2, ... ,77,2) . (6.8) 

The zij and thus the critical temperature depend only on the absolute effective charge, 
namely q = \ni — n2\/2. This property reflects the fact that the global U{1) sector of 
charge (77,1 + 772)/2 does not affect the thermodynamics in the large limit, since there 
are no degree of freedom which couples to it. We also observe that the critical equation 
is independent of the filling fractions sj and it is obtained by requiring the vanishing of 
the determinant 

(1 _ ^tot.(^\ _^tot.(„\ \ 
1 - 4- (x)) = - - = °' (^-9) 

where we have used that the matrix Vu is symmetric (^12*' — ^21*') ^^at z^ii' = 42*' — 
Zq"*' is the partition function in the trivial vacuum. This equation naturally splits into 
two simpler equations 

(a) : A_(x) = 1 - 4°*-(x) - z{f{x) = (6.10) 

(b) : X+{x) = l-z'°'-{x) + zlf{x) = 0. (6.11) 

The critical temperature is determined by the lowest zero of these two equations. Since 
— A_ = 2^12 > and A+(0) = A_(0) = 1, A_(x) reaches its zero at a smaller tempera- 
ture: in determining Xc we can then neglect A+(x). 

From the structure of the critical equation, A_(x) = 0, we can deduce two general 
properties of the critical temperature. First, the positivity of z\°2 also ensures that \-{x) < 
Xo{x) = (1 — Zq°^). This means that the critical temperature in a non-trivial monopole 
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background will always be smaller than the corresponding one in the trivial vacuum. 
Second, the function z*2*' decreases with the monopole charge q (in the interval x E [0, 1]): 
this implies that the critical temperature increases with the monopole charge. When 
q approaches infinity the value of the critical temperature becomes that of the trivial 
vacuum. Below we present a table for the critical temperature, where the behaviors 
described above are manifest 



q 


Xc 


T 


1/2 


0.085786 


0.407183 


1 


0.099771 


0.433863 


3/2 


0.103842 


0.441523 


2 


0.104567 


0.442884 


5/2 


0.104672 


0.443081 


3 


0.104686 


0.443107 


7/2 


0.104688 


0.443111 



Table 1: and Tc in the two sectors situation as a function of the relative monopole charge q. 

When the number k of sectors grows, the dependence of the critical temperature Tc on 
the relative monopole charges becomes quite intricate. However, some general behaviors 
can be anticipated. Consider, for example, a generic background of the form 



f = (rii, . . . ,ni,n2, . . . ,n2, 
where the induced relative monopole charges 



.rik, 



.rik 



(6.12) 



QiJ 



\nj - nj\ 



(6.13) 



are large, namely nj and nj are very different from each other. Then the Hagedorn tem- 
perature is dominated by the smallest charge and the off-diagonal terms associated to the 
other charges can be considered as small perturbations. The determinant is approximately 
given by 

(6.14) 



det{Vu) ^ {1 - zo)'~\{l 



^tot\2 
^0 ) 



- U, 



tot 



Exploiting what we have learned for the k = 2 system, the lowest transition temperature 
is an approximate solution of the equation 1 — Zq°^ — z''"^.^ = 0. 

Another interesting family of configurations is built by considering long sequences of 
sectors with equal length and monopole charge increasing by a fixed value q, namely 

f = (no, . . . ,no,no + q, . . . , no + q, no + 2q, . . . ,no + 2q, ,no + kq, . . .,no + kq). (6.15) 

When the number of sectors k goes to infinity, the Hagedorn temperature in these vacua 
approaches that of A/" = 4 on the Lens space S^/Zq in the sector described by a vanishing 
fiat-connection. For example for q = 1, a simple numerical analysis shows that Tc goes to 
that of pure A/" = 4, Tf =^ = -1/ log(7 - A^/3) ~ 0.379663 [9], (see table below). 



30 



K 








\J.\J00 I ou 


4071 84 


3 


0.079653 


0.395245 


10 


0.072873 


0.381820 


15 


0.072312 


0.380697 


20 


0.072098 


0.380267 


30 


0.071936 


0.379942 


60 


0.071833 


0.379736 



Table 2: Xc and Tc in the k sectors situation at 17 = 0. The vacua are labelled by 

h = diag(fc - l,...,k- 2, ...,0). 

Analytically, this result can be argued by noting that the matrix Vjj, of which we have to 
compute the determinant, is of Toeplitz type, namely a matrix in which each descending 
diagonal from left to right is constant. Consequently its entries do not depend on I 
and J separately, but only on the difference I — J. For this kind of matrices, when the 
dimension is large, the determinant is approximated by that of their circulant version |50j . 
This means that the smallest zero of the determinant can be found as a solution of 

oo 

1- J2 =0, (6.16) 

fc=— oo 

which is the smallest eigenvalue of the corresponding circulant matrix. In f l6.16p 2;*°* (/cq, x) 
is the single-particle partition function in the sector of charge kq. It is now possible to 
show that this infinite sum produces the single-particle partition function of the A/" = 4 
SYM theory in the trivial vacuum of S''^/Zq (see [39] for comparison). In other words (16.161) 
coincides with the critical equation for the TV = 4 SYM theory in the trivial vacuum of 

Finally we consider the addition of chemical potentials to a monopole configurations. 
Their introduction does not alter significantly the picture and a numerical analysis is 
given in fig. [2l 



6.1 Just above the critical temperature 

To understand what happens when we cross the critical temperature, we shall now focus 
our attention on the two-sectors configuration (16.81) . In this case, if we introduce the 
combination 

P± = ^(pi±P2) , (6.17) 

the action takes a diagonal form 

°° / 1 1 \ 

5 = 2 V -A_(x")p+„p+„ + -A+(x")p_„p_„, . (6.18) 

n=l ^ 

Above the critical temperature, A_(x) is negative and the dominant saddle-point is no 
longer realized by a flat distribution pin = p2n = (p+n = P-n = 0). In fact, as 
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Figure 2: Transition lines for three sectors vacuum: f = diag(i, 0, — z). Narrow lines corresponds 
to i ~ 1, thick lines to i = 10. The convention for continuous, dashed and dot-dashed are those of fig. [T] 
The qualitative behavior is the same for every number of sectors. 



the temperature is increased, the attractive term in the pairwise potential continues to 
increase in strength, so the eigenvalues become increasingly bunched together, occupying, 
at the end, only a finite interval / = [—6q,6q] on the circle (we arbitrarily choose the 
middle of this interval to be at 6' = for convenience). However, since A+(x) is still 
positive, we can safely assume that the new dominant saddle point satisfies 



P- 



2 (Pln ~ P2ri 



0, 



I.e. pi = p2 = P+ 



(6.19) 



In other words, the problem reduces to an effective one matrix model governed by the 
action 



s=2 d9de'p+{e)p+{e')J2 



n=l 



V. 



+ 2 / dede'p+{e)p+{e') log 



n 

e-e' 



cos{n{e - e')) 



sm 



(6.20) 



where the distribution function has support in the interval [— ^o? ^o]- In complete analogy 
with what we found in trivial vacuum case (16.211) . we have 



P{d) 



cos 



(I) 



vr sm 



'sm \ — 



9 

sm- , - 



with cos ( — 



X4x) - 



-. (6.21) 



Near the critical temperature, for T > Th we have the following expansion for the partition 
function 



F 

iV2 



-H 



A_(x) + 0((A_ 



+ 0{{T-THf), 



(6.22) 



T=T„ 
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which gives the characteristic first-order transition, aheady found in the four- dimensional 
model. 



7 Thermodynamics in non-trivial vacua II 

We discuss now the thermodynamical behavior arising when the second regularization 
scheme, considered for the fermions in section 4, is adopted. As previously derived, a 
non-trivial logarithmic deformation of the multi-matrix model (16. ip has to be considered 



k / oc ^ \ ^ 

X[[dU:]exV EE-^/°i-(^")TKf^7)Tr(f/?) X{det{Urr-'-^\ (7.1) 

1=1 \ IJ n=l / 1=1 

To illustrate the effect of the new interactions on the large dynamics, we shall make a 
very drastic assumption and we shall focus our attention just on the first winding, n = 1. 
With this choice the original matrix integral reduces to 

/ JjDf/,exp z\f{x)i:i{Ui)i:i{U\) JJdet(f/,)^"^-'3. (7.2) 
1=1 \ IJ / 1=1 

It is useful, as a first step, to introduce a set of k complex Lagrange multipliers A/ and 
the partition function can be written as 



^iilg-^ j \[dXjdXjeM-Y.^iNj\izj}{x)\j)x 
, " " (7.3) 

\\ I DUiexp (XiNiTi{Ui) + A/iV/Tr(f/j)) det{Ui)^'''~^ . 

T — 1 J 



Next we use the polar decomposition A/ = 7/e*°^ for each Lagrange multipliers. The 
phases e*"^ are then decoupled from the matrix integration by means of the change of 
variables Uj i— > Uje^"'' . This procedure yields the following integral 

(dettz,j)V J n "^^^"^"^ I - E NiNj^izY}{x)e^^^^-^'hj + i ^ Ni{Nni - Q)aj 

n / DUie^v [liNj{Ti{Ui) + Ti{U\))) det(t/,)^"^-«. 

r 1 ^ 



X 

7=1 



(7.4) 



In (I7.4p the group integrations over the unitary matrices Uj are completely decoupled. 
Each matrix integration corresponds to a Gross- Witten model [35] with a coupling 7/ 
and an additional logarithmic potential proportional to \og{det{U i)) . We remark that 
these kinds of deformations for unitary matrix models were widely considered in the 
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early eighties (see e.g. [271 [28]). The determinant operator was expected to act as an 
order parameter for the large N phase transitions characterizing this class of models 
[5T] : unfortunately, we cannot simply borrow the old results. In (17.41) in fact we have 
a new and decisive ingredient with respect to the original investigations: the power of 
the determinant is not a fixed number, but it grows linearly with A^. This last feature 
dramatically alters the usual large dynamics since the integral (17.41) is not dominated 
anymore by the same family of saddle-points of the familiar Gross- Witten model, as we 
will see in the following. 



7.1 Solution of unitary matrix model with logarithmic potential 

The phase structure of (17. 4p can be naturally studied along the lines proposed in |52j . 
We will first perform the integration over the unitary matrices and then the integration 
over the Lagrange multipliers. We will then start by studying the large properties of 
the reduced model 

Z(7,p) = j DUexp {'jN{Tt{U) + Tr([/^))) det{U)^P, (7.5) 

where Np is an integer, whose sign is irrelevant because we can transform Np into —Np 
by performing the change of variable U ^ . For this reason, from now on, we shall take 
p to be positive. The first important effect of the new logarithmic interaction concerns 
the small 7 behavior of ( 17. 5^ : differently from the Gross- Witten model (p = 0), where 
Z{'y,p) is finite as 7 approaches zero, here Z{'-f,p) vanishes as 7^ ^. This leading behavior 
is determined by expanding the exponential around 7 = and performing the integral 
term by term. The first non-vanishing contribution is fixed by the selection rule imposed 
by the U{1) factor present in U{N) and it is given by 



/ ]\T\N'''p f N-1 .. 

(7^6) 



where the constant C in the large limit is given by 

C = -^ ((log(4) - 3)p +{p+ 1)2 log(p + 1) - / log(p)) . (7.7) 

In other words, the free energy T{^,p) = log 2^(7, p) = N'^J-'o{'y,p) + .. of the present 
unitary matrix model starts, at leading A^^ order, with a logarithmic singularity similar 
to the one of the usual Penner model [531 El] 

J-o(7,p) = plog(27) + C + 0(7'). (7.8) 

This new behavior suggests that the usual strong-coupling expansion of (17. 5p might be 
radically different from that of the Gross- Witten model, which is simply given by e^'^'^^ . 
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To explore this idea, one could perform a full strong-coupling expansion and to resum 
the resulting series in the large limit; however the presence of the determinant factor 
much complicates this approach. Here, we shall choose a simpler path and consider a 
different expansion, peculiar of the present model, namely p very large. In this limit we 
can perform a semiclassical analysis on the integral fl7.5p : the relevant classical potential 
is, in this case. 



p i=i i=i J 



(7.9) 



The equations for the critical point are easily derived and solved (we will denote from 
now on 47^ = t) 

-ysin^i + z = 9, = isinh-^ (^-^ . (7.10) 

The semiclassical approximation is then obtained by expanding the classical action around 
the critical point up to the quadratic order 

(v^T+^-psinh-i _ N j^^/^j f^^^ + O (^7) , (7.11) 



with 



= -zsinh"^ ( ) . (7.12) 



We remark that this is a good approximation as long as \/jP~-\-t ^ 1: in this limit the 
gaussian integration covers the whole real line and the Haar measure over the unitary 
matrices becomes the usual measure over the hermitian matrices. We can easily perform 
the integration over the angles 6i and up to a constant independent of p we get 



Mt,P) = (^V^Tt -psinh-' (^-^^ - l/21og(v/p2 + t)^ 



p2 + t-plog| ^ + -^^ + 1 1 -l/21og(V7+t) 1 • 

For p large and t finite or small, we finally arrive to the following expansion 
M^) , .^ , t 1/1'' 



(7.13) 



Mt,p) =p { - log {p) - log(2) + Ij - 2 log(P) + - 



1 o /1\^ 1 o /1\^ „ / /l^ ^ 



32 Vp/ 8 \pJ Wp 



(7.14) 



This result is quite remarkable: we see that the above expansion reproduces exactly the 
large p limit of (17.81) and contains a infinite series of corrections in powers of t. Since 
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(17.141) holds also for small t, we must conclude that the strong-coupling expansion of 
our deformed Gross- Witten model leads to a non-trivial function of t and p, eventually 
encoding an intriguing modification of the p = result. 

It is quite easy to repeat the same analysis taking t large, exploring in this way the 
deformation of the weak-coupling phase of the familiar unitary model. In this limit we 
should obtain, at leading order in t, the very same result for the free energy as in the 
Gross- Witten case: again we could expect a non-trivial deformation due to the presence 
of p. Actually, performing the same steps as before, we get again (l7.13jM which expanded 
for large t gives 



^o{t,p) = - 




(7.15) 



We recognize in the first three terms the exact large result of the Gross- Witten weak- 
coupling phase: it does not come as a surprise, being the semiclassical approximation exact 
in this phase. As expected, we also observe an infinite series of corrections, depending on 
p, that modify non-trivially the usual spherical free energy of the weak-coupling phase. 

We do not expect, of course, that the above expansions yield the exact large A^ free 
energy: these results are semiclasssical, in the sense that we missed the contribution of 
the Vandermonde determinants associated to the measure over unitary matrices, that is 
essential in recovering the correct spherical free energy. Nevertheless they should capture 
the leading order behavior at large p or t of the complete large A^ answer, and also a 
certain series of subleading terms (as we will explicitly check in the following). 

These computations suggest an intriguing possibility: we observe non-trivial deforma- 
tions of both strong and weak-coupling expansion of the Gross- Witten model, involving 
complicated functions of p and t. It is quite natural to conjecture, at this point, that a 
unique non-trivial analytic function ^^{t^p) exists, reproducing for p 7^ both behaviors 
and being the large A^ free energy of the model. This is also suggested by the fact that 
the same free energy (17.131) describes smoothly either the large p or the large t region 
(see footnote [9]) . If this is the case, the presence of the logarithmic interaction would 
smooth out the third-order phase transition of the Gross- Witten model, the parameter p 
providing an analytic interpolation between the strong and the weak-coupling phase. 

In order to prove this idea, we have to solve exactly the large A^ dynamics: we shall 
exploit the beautiful relation between our model and the Painleve III system illustrated 
in [29] . In that paper the authors have shown that it is possible to construct an auxiliary 
function, 

a{t) = -t|log {{tNY'-^'l^e~'^'^'^Z{t,p)) , (7.16) 

^ The semiclassical computation really holds for \Jp'^ +t > > 1 and this condition is realized by taking 
either p or t large. Therefore we have to obtain the same free energy (17.131) in the large t case as well. 
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that satisfies, at finite N, the following non-linear differential equation 



- ^p^N^ + - 1) a'itfN^ + a\t) {Aa\t) - N^) {a{t) - ta\t)) + t^a^tf = 0. (7.17) 

In the large N limit, the spherical ansatz for the partition function Z(t,p) = e^^-^"*^*'^-' 
dictates the following scaling for the auxiliary a{t) 

a{t) = N^p{t). (7.18) 

Thus, at the leading order in A^^, we obtain a nice first-order differential equation for the 
reduced function p(t) 

- Atp'itf +{p' + t + Apit) - 1) p'{tr - pit)p'it) - = 0. (7.19) 

The analysis for small and large t given in (17. 8p and (17.151) provides two possible boundary 
conditions for the above equation: 

(s): p(t)|i=o = -|(p'+p); 
(w): p{t)\t^^ = {-^Vi. 

Since (17.191) is a first-order differential equation, these boundary values will correspond, 
in general, to two different solutions: the former, which satisfies (s), is denoted with 
Ps{t) and it is supposed to describe the strong-coupling regim^; the latter, Pw{t), obeys 
(w) and it is expected to hold in the weak-coupling regime. The two corresponding free 
energies J-'Q^{t,p) are then constructed by integrating the simple relation 



dJ^o'"it,p) fl p^ Ps,w{t) 



dt V 4 2t 



(7.20) 



which follows from (17.161) once we have used the spherical ansatz Z(t,p) = e^^'^o' 

The above simple picture works very well at p = 0, where our model reduces to 
the usual Gross- Witten model. In this case the differential equation becomes extremely 
tractable, factorizing into two simple first-order equations: the solution jFJ(t,0) and 
jF^(t,0) can be obtained explicitly and they exactly coincides with the well-known free 
energies of the model at strong and weak coupling. The condition J^Q{t,0) = J^Q{t,0) 
defines the correct critical value for the coupling constant {t^ = 1). When p ^ 0, the 
situation reserves some surprises as we shall illustrate below. 

As thoroughly described in appendix [Dl the general case can be solved exactly, in spite 
of the apparent difficult non-linearity of the differential equation. In particular there are 



^"in the matrix model language, 7 is conventionally identified with the inverse of the fundamental 
coupling constant. Thus small values oit = 47^ are in the strong-coupling region. 
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two relevant solutions, describing respectively the deformations of Ps{t) and pw{t) found 
in the Gross- Witten case. Integrating (17.201) we get a candidate J^^it^p) given by 

•^0 (i, P) = -\ (aog(4) - 3)p +{p + If \og{p + 1) 

while J-'Q{t,p) has the form 



2ptanh ^ (p + ^(j,~P^ P'w^ + 



(7.22) 

where the constant is given by 

/«, = + \p{{-3 + log(16))p - 2 log(p - 1) + 2 log(p + 1)). (7.23) 

Here p^(t) is the solution of the fourth order algebraic equation (]D.25p . which respects 
the large t behavior implied by the boundary condition (w). One can easily check that 
T^it^p) smoothly reduces, as p goes to zero, to the free energy of the Gross- Witten model 
in the weak-coupling phase, and accurately reproduces the semiclassical expansion ( 17.15^ . 
up to higher order terms in p2n^^n+m/2^ coming from the exact large solution encoded 
into the differential equation. It is also evident from (17.211) that J^Q{t,p) reproduces, in 
the limit of vanishing p, the Gross- Witten strong-coupling result. 

On the other hand, we already know that ^^{t^p), as given by (17.211) . cannot provide 
the right solution describing the small t regime! The large p expansion of (I7.2ip is quite 
boring and does not reproduce the non-trivial series (17.141) . obtained from the semiclassical 
approximation. On the other hand it is possible to show that jF^(t,p), for p 7^ 0, also 
satisfies the right boundary condition to describe the strong-coupling region (see appendix 
[D|) and, more importantly, correctly reproduces (I7.14p in the large p limit (up to higher 
order corrections in t^/p^""*"™, coming from the exact large solution of the model). 

We arrive therefore to the conclusion that the critical behavior of the standard unitary 
matrix model is completely modified by the addition of our logarithmic interaction. As 
long as p 7^ the system is always in a "weak-coupling" phase, described by the free 
energy ^(^(t,p): this solution has the correct boundary condition both at small and 
at large t and smoothly interpolates between them. We also identify J^Q{t,p) with an 
unphysical solution of the differential equation (I7.19P and therefore we neglect it. The 
situation drastically changes for p = 0: it is possible to show that, starting from J^Q[t,p), 
the limiting behavior changes discontinuously at t = 1. On the other hand, taking p = 
at level of the differential equation (17.190 . the strong-coupling phase is instead encoded 
into the solution Ps(t). 



t p 
-p^\og{p)) + j^^-—^--\ogit) , 

(7.21) 
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7.2 Phase-structure in non-trivial vacua 

In this subsection, we shall explore the consequences of the previous results on the phase 
structure of the theory. After having performed the integration over the unitary matrices 
in the deformed Gross- Witten models, we are left with the integration over the Lagrange 
multipliers 

„ k / k \ 

/ J] d^idai exp Y.^tsi{ni - q)aj + s]To{-ij,p)) -N^Yl s/Sj7/^7i(a;)e*^°^"°^Vj 
j=i \ 1=1 IJ / 

(7.24) 

with q = Q/N. Since N is large, we can perform this integral in the saddle-point ap- 
proximation as well. The saddle-points which dominate this integration are determined 
by 



2 ^ sniZjj'jj sin(aj - a/) + i{nj - g) = 



7=1 
k 



(7.25) 



-2 ^ s^iZj} cos(aj - ai) + sjJ^q{'Jj,p) = 0. 



1=1 



To be concrete, we shall consider only the case k = 2: here the relevant combinations of 
the parameters are given by ni — g = 52(^1 — n2) = and n2 — q = —Si{ni —n2) = —Sin 
(ra = rii — 77-2 > 0). The first equation in (17.251) then produces two conditions 

272-22i"'^7i sin(ai — 02) + in = and 2'yiz^2l2 sin(a2 — tti) — in = 0. (7.26) 

These two equations are obviously equivalent and they are solved by 



n , . . I n^ 



sin(ai - 02) = -i- r]— =^ cos(ai - ^2) = ±W 1 + 77 n — r^- (7.27) 

272^21 7i V 4(72^21 7i) 

Substituting this result into (17.251) . the second equation provides two relations, which 
determines 71 , 72 



2si7?2i/ T '2s2\l (72^2/71)^ + X + ^171-^0(71, P) = 0, 

^ (7.28) 



T 2Siy (72^21^1)^ ^ T ~ ^^272^11^ + S272^o(72,P) = 0. 

In the following, we shall further simplify our example and we shall choose two sectors of 
equal length, namely we shall set Si = S2 = 1/2. Then, by taking the difference of the two 
equations and using the fact that J^o{t,p) is a monotonic function, one can immediately 
show that 7i = 72. We remain with just one equation, which determines ti = 47^ 



^{z^2'hy + n^- l^iY + 2ti^^(ti,p) = 0, (7.29) 
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which is conveniently rewritten in terms of Pw{ti) as follows 



= ±y ^{z^iti)' + n' + j{l~ ^1') +n = 2 (p^(ti) + ^(n + 1)) . (7.30) 

When ti runs from zero to infinity, the r.h.s of fl7.30p spans the same region. Thus a 
necessary condition for having a non-trivial solution is that the l.h.s. of fl7.30p is not 
negative definite. Let us discuss the first equation: has the following properties 

/+(0) = 2n, /;(0) 

= ^ f 

Moreover we have that for large t 

/H-(t)^^ ""+f + (7.32) 

^ Zii -\- Z\2 

We immediately conclude that for temperatures near zero (x <^ 1), /+(t) is always de- 
creasing: for T < Tf/, where Tu is the Hagedorn temperature defined by the equation 

^11 + ^12 = 1 (7.33) 

as in fl6.10p . we have that det(l — -z) > 0, implying that /+(t) never vanishes for t > 0. 
Therefore in this range of temperature there is always one solution to the saddle-point 
equation at t 7^ 0. At the Hagedorn temperature Tu the function /+(t) is still decreasing 
but becomes positive definite, asymptotically approaching the value n. Above we see 
that /+(t) develops a minimum at finite t and then becomes monotonically increasing. 
The minimum disappears at the temperature T2 defined by 

det z = zn (7.34) 

and the function becomes monotonically increasing for any T > T2. 

In spite of these changes of behavior, one can check that there is always one solution 
to the saddle-point equation as shown, in different regimes, in figure [3l Moreover the 
position of this saddle-point changes smoothly as function of the temperature (fig. HI). 

Let us examine the second saddle-point equation, the one involving f-(t). We can 
repeat the same analysis: the main conclusion is that for temperature T < T2 we see /_ (t) 
being monotonically decreasing and therefore, because /-(O) = 0, there is no solution for 
t 7^ to the saddle-point equation. We notice that t = is not acceptable because of 
f l7.26p . For T > T2 it is not easy to see analytically if f~{t) provides new solutions to 
the saddle-point equation: we have done a numerical study, showing that a new solution 
appears for x > 0.212352. However, the resulting free energy is always subdominant with 
respect to the other one as illustrated in fig. [5l So the solution associated to /+ is the 
only relevant saddle-point in the large limit. 



1 det z 
2 



-An' 



det z 
(det z 



Zll) 



det z 



^12 



det(l 



(7.31) 
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Figure 3: Plot of /+(t) for different values of T and n = 1, q = 1/2. Going bottom-up, the solid lines 
illustrate the behavior for T <Th,T = Th (lower thick line), Th < T < T2, T = T2 (upper thick line), 
T > T2. The dashed line is the r.h.s. of (|7.30p as a function of t. 



t 




Figure 4: Saddle point as a smooth hmction of the temperature x for n = 1, q = 1/2. The graph 
covers all the different regimes, the Hagedorn temperature being xh = 0.0857864 and T2 corresponding 
to X2 = 0.115493. 
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Figure 5: On the left side the saddle-pomts (for n = 1, q = 1/2) m terms of p'{x) associated to /+ 
(continuous hne) and /_ (dashed line) are shown. At x ~ 0.212352 the solution intersects with the 
t = unphysical solution (dotted line). On the right the free energies for both cases. 

We conclude therefore that within our approximation, that consisted in taking just 
the first winding in the matrix model action [n = 1), we have always a non-trivial saddle- 
point giving a free energy Fb = log Zb of order A^^. Moreover this saddle-point varies 
continuously with the temperature: in particular at the Hagedorn temperature Th-, repre- 
senting the point of the first-order phase transition in our first regularization scheme, the 
free energy remains smooth and no discontinuous behavior appears in this second scheme. 



8 Conclusions and future directions 

In this paper we have studied the maximal supersymmetric gauge theory on R x S*^, 
with particular attention to its thermodynamical properties in the limit of zero 't Hooft 
coupling. In the case of the trivial vacuum, we found a behavior similar to the parent 
four-dimensional theory, with a first-order Hagedorn transition separating a "confining" 
phase from a "deconfined" one, with non-trivial expectation value for the Polyakov loop. 
We have repeated the analysis for monopole vacua and we have apparently different 
behaviors, depending on the regularization procedure: this actually reflects the particular 
choice of the fermionic three-dimensional vacuum, that is related to generation of Chern- 
Simons terms when monopole are present on the sphere. We have presented two opposite 
choices, both allowed at quantum field theory level, generating different unitary multi- 
matrix models describing the thermal partition function. The critical behaviors we found, 
under suitable assumptions on the relevant contributions at small temperature, are very 
different: in particular we have observed that no Hagedorn transition seems to be present 
within our second regularization choice. Further studies are surely necessary to elucidate 
the situation: first of all we expect that supersymmetry should play a role in order to 
distinguish between the different regularization choices and consistency with the SUSY 
algebra could probably select a preferred "vacuum charge" . On the other hand the relation 
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with the gravitational duals should also be investigated to provide a physical interpretation 
of the Casimir energies and of the Chern-Simons contributions. Apart from solving the 
puzzles arisen in this paper, there are a lot of potential interesting developments involving 
the study of the A/" = 8 three-dimensional supersymmetric theory considered here. It 
would be important of course to determine the nature of the phase transition beyond 
zero 't Hooft coupling and to discuss the issue of exact decoupling limit using chemical 
potentials, in the spirit of [20l [23| [2^. We also plan to consider the phase diagram in the 
presence of background scalars as in [551 EE] . More generally one could try to explore if 
some remnant of four- dimensional integrability persists in three dimensions and to make 
some quantitative connection, in the strong-coupling limit, between the gauge theory and 
its gravity dual. It would also be very interesting to study BPS Wilson loops on S*^: in four 
dimensions there have been exact results for particular classes of loops, the computations 
reducing to matrix integrals [571 EH] • It is natural to ask if a similar phenomenon takes 
place in three dimensions too. 
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A Conventions and supersymmetry variations 

Before discussing in more details the supersymmetry variations considered in section[2l 
we shall briefly summarize our conventions and identities on F-matrices. 

Metric and gauge conventions: The metric is taken diagonal and with Minkowskian 
signature: TjMN = { — , +, +}• The capital letters M,N, . . . will span the ten dimen- 
sional spacetime indices (0, 1, . . . , 9), while the Greek letters /i, z/ . . . will denote the three 
dimensional spacetime indices (0, 1, 2). The indices i,j, k are associated to the directions 
(1,2) along the sphere S*^, while the directions (3, . . . , 9) transverse to 5*^ are indicated 
with m, n, . . . . Finally a special index notation is also reserved to the set of directions 
(4, . . . , 9) for which we shall use the overlined letters rh,n, . . . . 

The gauge fields A = A°'t°' are taken to be hermitian and the generator are nor- 
malized so that Tr(t"t^) = The covariant derivatives are then defined as follows 
= — ig[Afj_,-], where is the geometrical covariant derivative. In general we 
shall omit the trace over the gauge generators in our expressions, unless it is source of 
confusion. 
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Some useful F-identities: For convenience, here we have collected some F-identities, 
which are useful in checking the super symmetry invariance of the Lagrangian of our model: 







-pO-pjfc-pO _ 


_-pjk 




pOpOjpO ._ 


= FO^ 


pjpjmpi 


0, 








■pOpOmpO 


pOm 




12'pmn 




pmn 









(A.l) 

Summation over repeated index is understood. Here F*^ denotes the ten dimensional 
matrices, while the symbol is used for the three dimensional Dirac matrices. The 
symbol r^^i*'^^ defines the completely antisymmetrized product of the matrices F^^, 

Three-dimensional fields: The scalar field <^*-' is antisymmetric in i, j, which are 
SU{4)ji indices and it satisfies reality condition: 

<P'' = i<P.,y = le^'''<Pki. (A.2) 
It is defined in terms of the old fields 0^ by the relations: 



h4 -t- (P23 , -(Pl3 -t- 924 , 912 "t" 934 



. .014-023 , .013 + 024 , ."012 + 034 



(A^3) 



The spinor fields Aj (again, i is an SU{A)ii index) denote the Dirac spinors in D = 3 orig- 
inating from the dimensional reduction of ipM, while describes the three-dimensional 
gauge field. 



A.l Super symmetry variations 

In this appendix, for completeness, we shall write the conditions for the vanishing of the 
variation at the order a and at the order a^. At the linear order the complete variation 
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can be summarized by the following table: 



Term 


Condition 


2Re{a^7[0M,0n]V^r™"ri23e} 


B + 2 + P + M = 




2Ee{a(?[03,</>nr]^r3-ri23e} 


2i3 + 4 + 2P + G - 2M = 





2Re{mDo03?rO3ri23e} 


4 - 2i3 + P + G - 2M = 





2Re{aiL>o(/>MV'r°'^ri23e} 


A — 2B + P + 2M = 




2Ee{aiA03?r*3ri23e} 


2S + P + G + 2M + A^ = 





2Re{mA0m?r*™ri23e} 


2i3 + P - 2M = 




2Re{aiFoi^r°Ti23e} 


-2B -2M = 






i3-2 + M+ f = 





(A.4) 



There are eight different kind of terms, listed in the first column, and they must vanish 
separately: this leads to the conditions in the second column. 

At the quadratic order in a we have simply 



Term 


Condition 


2Re{ia20_^r^^} 


-2V+ (2- ^)P + MP = 


2Re{ia^(^3i^T^^} 


-2{V + 1^) + (2 - £)(P + G) - M(P + G) = 



(A.5) 



B Computing the one loop partition function 

Here we give all the details of the calculation of the partition function in a monopole 
background. For the free model the one-loop contribution of each field is a functional 
determinant, giving the single-particle partition function. 



B.l Computing determinants: the master-formula 

We illustrate our regularization scheme: readers who are not interested in these details 
can take (IB.SP and (]B.17p as main results, and skip to next subsection. 
All the determinants appearing in the evaluation of the free partition function contains, 
as a key ingredient, the evaluation of the following infinite product 

2j+P 



oo oo 

n n 

j=0 n=— oo 



(j + v) + 



47r^ 



In + w) 



(B.l) 
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This quantity is divergent and it must be regularized. Here, we shall adopt the standard 
C— function regularization and we shall define 

E(r?,p,/3,^):=e-f'(°), (B.2) 

where 

2j + P 



oo oo 



Notice that fIB.SP defines the function (^(s) only for |s| > 1. In order to compute C'(0), we 
have to consider its analytical continuation to a neighborhood of the origin in the s-plane. 
This is achieved through a standard technique: firstly, we shall use the Mellin-Barnes 
representation and subsequently we shall perform a Poisson-resummation in n 

100 OO ft r, 



^^4— 5^(2j + p) /rft f-|.-*0■^-'')%"^-2-- 



j=0 



2 e " " e « 



/3r(.-|)^ (2j + p) , 



/3r(.-|: 
2v^r(s) 



(2C(2s-2,r/)-(2r/-p)C(2s-l,r/)) 



The only contribution to C'(0) ( 1B.4I) arises when the derivative acts on l/r(s) since 
this quantity vanishes as s approaches 0. We obtain 



oo oo 



. , , , , v^v^ 2e-"^(^'+'')(2? + p)cos(2ra7ru;) , , 
C'(0) = -/3(2C(-2, v) + {p- 2r/)C(-l, ^)) + E E -^—^ ^ -■ 



n 

j=0 n=l 



From the final expression ( IB. 51) we can deduce two equivalent representations of this result, 
which are both useful for our goals. Firstly we can perform the sum over j, which yields 

C'(0) = (IbM + hp- 2r))B,{r))) + "/"^^^^^ cos(2n™), (B.6) 

with X := and Bk{ri) being the Bernoulli polynomial. Next, we shall define the 
"single-particle" partition function 

(1 -x) 
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and finally write 

log (S(r/, p, P,w)) = -P C^Bsiv) + lip- 2v)B2iv)) - 2 """^^^"""^ cos{2nnw). 

^ ' n=l 

(B.8) 

This representation will be the most natural when discussing the matrix model and the 
position of the Hagedorn transition. 

Alternatively, in (IB.Sp we can first sum over n 

2 1 °° 

C'(0) = f3{-B,ir]) + -(p - 2v)B2iv)) - 5^(2j + p) (log (l - zx^) + log (l - zx^)) (B.9) 

j=0 

where z := q-Pv+'^^-^w^ jf define 

CO oo 

r/(z,x) := and M{z,x) ■.= Y[{l - zx^ , (B.IO) 

i=o j=o 

we can recast the above result in a very compact form 

S(r/, p, /5, tf;) = e-/3(|^3 W+i(P-2^)B2(»,)) |^(^^ |2p |4 _ ^g^^^^ 

This second representation will be the most suitable when discussing the high temperature 
behavior. In this limit the leading contribution is encoded in the function Fp{z,x) 



oo 



Fp{z, x) = Y,{2j + p) log (1 - zx') . (B.12) 

j=0 



The X 1 behavior is transparent by rewriting Fp{ 

°° z^ 

Fp{z,x) = -Y,— 

' ^ m 



m 

m.=l 



^ 1 2 

p-2 + 



(B.13) 



and expanding in /3, at fixed z, we get 

2t-/\ Pt-/\ /P~2 5 



= -77^ Li3(--) - ^ Li2(--) + + log(l - z) + 0(/3). (B.14) 



To recover fl5.13p . where the contribution of chemical potentials to the high-temperature 
limit has been presented, we need further expand Lis (2;) for z ^ 1: we are interested in 
the case when w = and w = 1/2, appearing respectively in the bosonic and fermionic 
case, and with zero flat-connection {z = e~^, y ^ 0) 



U,{e-y) = C(3) - !^ y + Q - 1 log 1/2^ y' + 0[y') 
Li3(-e-^) = -^C(3) + ^y-\ log(4) + 0{y% 



(B.15) 
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Fermionic zero modes: In order to compute the contribution of the fermion zero 



'2ti . 

— in + w) 



If we adopt the zeta 



modes, we need to compute the product = Y\ 

n=—c 

function regularization as before, we are led to compute the following accessory sum 



G{s) 



(27r) 



.p(C(s,w) + e'''X{s,l-w)). 



Then 



log(i?) = -G'(0) 



n=l 



n 



(B.16) 



(B.17) 



B.2 The scalar determinant 



Let us discuss the solution of the eigenvalue problem (14.61) . Since our background is 
static, we can factor out the time-dependence in the eigenfunction by posing 9, 0) ~ 



Then the eigenvalue problem in the Weyl basis (14. 7p takes the form 



aSroots L 

N-1 



47r2 / 



P 

A0Q„ + —(pan + P^q^ 



ceYcen 



i=l 



47r^n^ /i~ 

-^<P^n + —(pin 



Min 1 H 



^ 1=1 aGroots 



(B.18) 



where A denotes the geometrical Laplacian for a scalar on the sphere. The symbol 
A instead represents the geometrical Laplacian in the background of a U{1) magnetic 
monopole of charge g^- This Laplacian is constructed with the covariant derivative 



D„ = V, 



iqaA 



(B.19) 



where is the geometrical covariant derivative. In (IB.ISP the components along the 
different directions in the Lie algebra do not interfere and we can consider them as inde- 
pendent. This allows us to split the original eigenvalue problem into two subfamilies, we 
have: (a) N{N — 1) independent eigenvalues coming from each direction along the ladder 
generator 

' A0a„ + ^(pan + P^q^fpan = Kn(pan, (B.20) 



n 



2tx 



and (b) — 1 independent eigenvalues coming from the directions along the Cartan 
subalgebra 



An + -T^m - A0i„ = XinCpi, 



f3^ 



(B.21) 



To begin with, we shall focus our attention on the family (a), since the family (b) can 
be obtained from (a) as a limiting case for a^, qa ^ 0. The solution of the eigenvalue 
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equation flB.20|) can be translated into an algebraic problem if we introduce the angular 
momentum operator in the presence of a U{1) monopole of charge Qa- Its form [59] is 



(B.22) 



Here are the Cartesian coordinates of a fiat R'^ where our sphere S"^ is embedded. In 
terms of this auxiliary operator, the kinetic operator in (1B.20P takes the form 



47r2 

/52 



n + 



4 



(B.23) 



Thus our task is reduced to finding the eigenvalues and the eigenf unctions of this dressed 
angular momentum operator (L*^"))^. Its spectrum!^ was determined thirty years ago by 
Wu and Yang [SI] and it is formally equal to that of the usual angular momentum: the 
eigenvalues are ja{ja + 1) and their degeneracy is 2ja + 1. What changes is the range 
spanned by the index ja, which now is Ig^l + 1, 1^^1+2, ■ ■ ■ . Putting everything 
together the spectrum of the kinetic operator flB.201) turns out to be 



A. 



47r2 



n 



with = \qa\ + l, |ga|+2 



(B.24) 



and each eigenvalue has degeneracy 2ja + 1. Notice that the spectrum does not depend 
on the sign of g^. The contribution of the family (a) to the effective action is given by 
the infinite product 



r*(„) = log 



n n n 

a&oots j^=\q^\ n=— oo 



An' f (3a, 



2ti 



2-1 2jc+l'^ 



(B.25) 



which is easily computed by using the results of appendix IB. II (with p = 1 + 2\qa\, 

we obtain (14. 9 p and (I4.10p The contribution 



7] = 1/2 + \qa\, w = ^^)- Setting x = 
of the family (b) is then obtained from the above results by setting Qo 



0. 



B.3 The vector/scalar determinant 

The eigenvalue problem for the coupled system (03, A) can be simplified by choosing the 
gauge-fixing (14. lip . This choice allows us to cancel some of the mixed terms {(psA) in the 
Euclidean quadratic Lagrangian and to obtain 



,(2) 



A.D^D^A-' + R^.^A'^A-' -iF,,[A'',A'^] - [A„ 03][^^ 03] + 



2 j,2 



2-^(f),e''''%D,A,. 



[B.26) 



^^The eigenfunctions are also known and they are given by the so-called monopole harmonics Yqjm {0, (p). 
They are a straightforward generalization of the usual spherical harmonics, but we shall not need their 
explicit form here. We refer the reader to [59] for more details. 
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Then, the following coupled eigenvalue problem 

- 003 + /iVs + [03, [03, 03]] - /iv^ep.AA;^^'^^^ = A03, (B.27) 
-UA, + R^,A^ + z[F,^, A'^] + [03, [03, A,]] + /iV^epA.A;^D^03 = \A,. (B.28) 

Since both the geometrical and the gauge background are static, the time-component of 
the vector field uj = k^Ap = Aq decouples completely from the above system. It satisfies 
the massless version of the scalar equation studied in IB.21 namely the eigenvalue problem 
associated to this component is 

-□^0^+ [03,[03,^]] = Aoo;. (B.29) 

We shall forget about u since its contribution is cancelled by the ghost determinant. We 
are left with the system given by flB.27|) and (IB.28P where the indices run only over space. 
We expand the coupled system (]B.27p and ( ]B.28|) in the Weyl basis and we factor out 

27vin ^ 

the time-dependence of the eigenfunctions: 6, 0) ~ Ann{6, 0) e and 03(t, 6, 0) ~ 



■>3n 



{6, 0) e . Along the directions associated to the ladder operators Ea we find 

— AAjan -|- m^Aian + '^f^'^ Qay/g^ijA-j^j^ + fJ^'^QaAian + l^^fd^jiD-^ 4''ian = ^anA-iani , , 

^ . (B.oUj 

where = (2TTn / (3 + a^)'^ + n"^ . In the first equation, the symbol A denotes the Laplacian 
on vectors in the background of a monopole of charge g^, while in the second represents 
the Laplacian on scalars. Along the Cartan directions we shall again get the system (1B.30P 
but for Qa = 0. 

To find explicitly the spectrum of system (lB.30p . it is convenient to decompose our 
vector Aian in its self dual part Af^^ and anti-self dual part v4~,„. Consequently we shall 
introduce the differential operators mapping ( ant i-) self dual vectors into scalars and 
their adjoints, mapping scalars into ( ant i-) self dual vectors. They are defined by 

(B.31) 

where D as in (IB. 301) stands for the covariant derivative in the background of a monopole 
of charge g^. In terms of these operators, the system ( IB. 301) takes the form 



2 2 f2 \ 

„2 2 f)2 \ 

+ + - |oi-»03«n = (B.32) 



■'San • 



Here we have dropped the index i because it is immaterial and we have set = 1^^ + jj?. 
At first sight the eigenvalue problem might appear cumbersome, but in this representation 
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it is quite simple to provide a basis where our problem reduces to diagonalizing an infinite 
set of three by three matrices. In fact, let us take qa > 10 and consider the following 
basis for scalars, selfdual and anti-selfdual vectors on the sphere 

= of ^V,^,^ for J > g„ + 1, e"^ = 0^^^%^jm and eJ^J, = Y.^.^n for j > g„. 

(B.33) 

Here Yg^j^ are the monopole harmonics, namely the eigenf unctions of the angular mo- 
mentum (lB.22p . For the anti-selfdual vector we have to add also 2{qa — 1) + 1 elements 



coming from the zero modes of O^^^of"''. We shall denote them as _x)m- -^^^ ^ 
tailed proof that ( 1B.33I) with the addition of the zero modes is a basis, we refer the reader 
to [60], where the following two useful identities are also shown to hold: 



OtlPtl)^% = -qIt qa)Y'^-^"' = + l)-qlT q.))e% (B.34) 



and 



0g/0geij„, = 0g/0g0g/ei% = + 1) - T g«))ef^. (B.35) 



Because of ( 1B.34P and ( IB. 351) and the definitions (lB.33p . e^-^ and for fixed j > go- + 1 



and fixed m generate an invariant three-dimensional linear subspace for the eigenvalue 
problem (]B.32p . The original problem can be then separately solved in each subspace. 



where it reduces to diagonalizing the following three by three matrix 

^ml - + j{j + - -n^ 

ml-fi^+j{j + l)fi^ + q^fi^ -fi^ I . (B.36) 

-4 (-9a - 9" + J(i + 1)) {-qI + Qa + j{j + 1)) ml + j{j + 

The three distinct eigenvalues of this matrix are given by 
A+ = 4+jV, A_=£2 + (^- + i)V, X3 = il + j{j + l)fi^ withj >g, + l. (B.37) 

For j = Qa self-dual vectors do not exist and the invariant subspace is generated only 
by e~"^ and ej?"„. Instead of ( IB. 361) . we have the two by two matrix that is obtained 



from ( IB. 361) by dropping the first row and the first column. Its diagonalization produces 
the following two eigenvalues A„ = £^ + (g^ + 1)^ and = + QaiQa + !)• Finally, we 
have to consider j = — 1- In this case, we are left with a one-dimensional invariant 
subspace generated by e^^_-^^^. The eigenvalue is simply A„ = + g^- Summarizing we 
have = i1 + (j -|- for j > — 1 and X3 = il + j{j + for j > q^ so that we 

have extended the range of existence of the eigenvalues (]B.37p . The degeneracy is always 
2j + l. 

In the following we shall neglect the family with eigenvalue A3, since its contribution 
is cancelled by the ghosts. We shall just consider the first two families X±, which instead 



-•^^The case < — 1 is obtained by exchanging the role of self-dual vectors with that of anti-self dual 
vectors. The value qa = ±1/2 and qa = will be discussed separately. 
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yield the actual vector determinant. The contribution of A+ is obtained from the results 
of appendix IB. II by setting w = [3aa/ (2vr), rj = qa + 1 and p = 2qa + 3 



aSroots 



n=l 



n 



with 



(3 - x) 2go 



fB.38) 



(B.39) 



1^1 — x)^ 1 — X 

The contribution of A_ is instead obtained setting w = (3aa/ (2vr), i] = and p = 2qa — 1: 



QGrootsX 



18ql + 10q^-l)-2j2 



^vect. 



X 



Jnf3aa 



n=l 



n 



with 



vect. 



x(l + x) 



2qa 



X 



(B.40) 



fB.411 



[l-x)' 1 

When adding these two contributions, we obtain fl4.12p and (14.131) . For what concerns 
the non-negative values of the monopole charge, there are still two cases to be considered: 
qa = 1/2 and = 0. In both cases, the elements of the basis coming from the additional 
zero modes of the operator Q^^^^O^^^ disappear [60]. For = 0, in the basis (1B.33P the 
element ej" with j = = is absent. The net effect is to reduce the range of the 



existence of the eigenvalues A_ 



+ U + 1)V' to 3 > q, 



eigenvalues Ag 



+ j(j + l)p^ to j > 1 for qa 



for qa = 0, 1/2 and of the 
0. By recomputing the contribution of 



A_, for = 1/2 we get the same results flB.40p and flB.4ip when we use the appropriate 
values for p and r] in appendix IB. II r] = 3/2 and p = 2. This does not happen, instead, 
for go, = : by using r] = 1 and p = 1, we get (14.141) . 



B.4 The spinor determinant 

In determining the contribution to the total partition function of the spinors Aj, we shall 
follow closely the steps of the previous appendix. The fermion kinetic operator expanded 
around the background (13. 3p has the following eigenvalue problem 

- t^^D^X + ^[03, A] + ^^7°A = pA, (B.42) 

where we dropped the S'f/(4)^ index since all the components give the same contribution. 
Expanding the matrix- valued field A in the Weyl basis and separating the time-dependence 
we get 

A = ^inH' + E ^-^°) e-¥("+i)*. (B.43) 

\ £=1 QfSroots / 

The only real difference with the scalar and vector cases is that fermions have antiperiodic 
boundary conditions along the time circle. The usual procedure will, in turn, disentangle 
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the different components along the Lie algebra and it will divide the eigenvalue problem 



(IB.42P into two subfamilies. As in the scalar case, we have: (a) A^(A^ — 1) independent 
eigenvalues coming from each direction along the ladder generator 



-7 



-^(n+-|+a. 



4 



A. 



iYDiXan + iMaKn = PanKn, (B.44) 



and (b) — 1 independent eigenvalues coming from the directions along the Cartan 
subalgebra 



en 



-7 



27r / 1\ u' 
— \ n + - ] - i- 



in- 



(B.45) 



In f lB.45P the symbol V denotes the geometrical covariant derivative on spinors while Di 
in (]B.44p is the covariant derivative in the background of a U{1) magnetic monopole of 
charge g^, i.e. 

A = di + ^rf S,fe - tq^A. (B.46) 



We shall first consider the family (a). The problem of diagonalizing the operator ( ]B.44[) 
can be solved algebraically by exploiting the unitary transformation U = e^^^'^ea'^"'^ 
fact, after performing this transformation, the operator (1B.44P becomes directly related 
to the total angular momentum J^"^ = L^"-* + | in the monopole background 



27r 
J 



n 



4 



(B.47) 



Here f stands for the usual radial unit vector in three dimensions while cTj are the Pauli 
matrices. In flB.47p . the operator S is the sum of three contributions. There is a reduced 
Dirac operator 

25(a) ^ ^e^jfcf .(x^. =i^f..a + /ie*^'=f,a^Lf ^ , (B.48) 



which is the standard two-dimensional massless Dirac operator in the presence of a 
monopole, but written in an unusual basis. Then we have a "chiral" mass term pro- 
portional to (a ■ f), which plays the role of the two-dimensional 75 (we have in fact 
{{a ■ f), !l)''"^}=0). Finally there is a constant shift proportional to the charge qa- 

Now, we can focus our investigation just on the operator flB.48p . since the spectrum of 
flB.47p follows from that of For each eigenfunction il) of S'-"^ with eigenvalue p 7^ 

there exists another eigenfunction (a ■ r)il) with eigenvalue —p. The possible values of p 
can then be computed by considering the eigenvalues of (S^""-*)^. This operator has the 
following simple form 

1 



p 



{J' 



a)^2 



+ 4 



(B.49) 



and it is diagonal on the basis ipjm± of the total momentum eigenfunctions, which satisfy 
(cr ■ f)ipjm± = il^jm^:- The cigeuvalues are p^^, = p^((j + 1/2)^ — g^). The positivity of the 
operator (S)*^"))^ imposes (j -|- 1/2)^ — g^ — 0; ^^"^ ivxn j > \qa\ — \. The degeneracy 
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of each eigenvalue is 2(2j + 1). On this basis, the operator D^"^ is also diagonal and it 
possesses the following spectrum 

'S)^°'^ljjjm+ = Pja1pjm+ and ^^"''^ijjjm- = -pjalpjm-- (B.50) 

In ( 1B.50P each eigenvalue has degeneracy (2j + 1). The above analysis does not directly 
extend to the kernel of the operator D^°'\ which is obtained for J = |q'a| — |- These zero- 
modes can be classified by using the eigenvalues of the operator {a ■ f): we shall denote 
z/± the number of zero modes with eigenvalue ±1. Then a simple application of the index 
theorem shows that = \qa\ — Qa and z/_ = \qa\+ Qa, namely for positive qa we have 
only zero modes with negative chirality and viceversa. 

We now turn back to the problem of diagonalizing the operator S defined in flB.47p . 
The operator S on the basis provided by the eigenvectors of D^"^ is not diagonal. However, 
on the subspace spanned by the eigenfunctions of non-vanishing eigenvalue, it factorizes 
in an infinite series of two by two matrices. Each matrix acts on the space generated by 
the eigenfunctions ipjm± and it has the form 



^ J- 1\ - n ^nH 



2-K 



n 



13 V ' 2 



(B.51) 



Since we are only interested in the determinant of the operator S, we shall not really need 
to convert this matrix into a diagonal form, but it is sufficient the evaluate its determinant 

/i'(j„ + l/2)2 + — r2 + - + ^-z^ withj, = |g,| + -,|g„| + -,..., (B.52) 



/32 V ■ 2 2n Stt ; 2"""' 2' 

and to recall that there are 2j + 1 determinant with the same value. Then by using the 
master formula of appendix lB.il (with p = 2 + 2\qa\, rj = 1 + \qa\ and w = ^ + — iff), 
the contribution of this part of the spectrum gives fl4.17p and 04.181) . On the kernel of 
D^°'\ the operator S is instead diagonal and it has the following spectrum 



_ 

Pna+ - J 

_ 271 

Pna- - J 



1 Ptta , .Ppqa 

-n — - — h I- h i- 



2 27r Svr 2n 

1 (3aa .Pp , .Ppqa 
n + - + — h 2- 



with degeneracy \qa\ + qo 



with degeneracy \qa\ — qo 



(B.53) 



2 27r Svr 27r 

Now we have to deal with the regularization ambiguity discussed in section 14. 4[ In our 
case, all the different choices for the cuts in the s-plane are encoded in the two following 
situations: 

(I) we regularize the determinants associated to the zero-modes" of negative and 
positive chirality by choosing opposite cuts in defining the complex power (one on the 
real positive axis and the other on the real negative axis). With the help of appendix 
f lRTll . we then obtain fICTD : 

(II) we regularize the determinants associated to the zero-modes" of negative and 
positive chirality by choosing the same cut in defining the complex power. A similar 
analysis yields fl4.20p . 

The appearance of a Chern-Simons term for case II and the total fermionic contribution 
to the effective action for both cases are discussed in section 14.41 as well. 
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C U(l) truncation of A/^ = 4 super Yang Mills 



In this appendix we show that the previous results can be easily recovered from A/" = 4 
super Yang Mills theory on IRx S*^ by a suitable U{1) projection which gives the maximally 
supersymmetric theory on R x 5^. 

The single-particle partition function in the representation R, z^{x), is given by 



where E is the energy eigenvalue subtracted of the Casimir energy, which can be derived 
for example with the procedure described in the body of the paper. The eigenvalue E can 
be computed most directly by noting that the Laplacian on the sphere may be written in 
terms of angular momentum generators which can be diagonalized by means of generalized 
spherical harmonics on S^. The isometry group of is 50(4) ~ SU{2)i x SU {2)2 and we 
will need the spherical harmonics for scalars, vectors and fermions, which will be denoted 
by Sj^rn,fh{P), V,-,m,m(^) and Fj^m,jn(f^), respectively. Here m and fh are the eigenvalues 
of .P and for SU{2)i and SU{2)2 and VL represents the coordinates of S^. We follow 
here the notation of [30]. Having determined the single-particle partition functions on 
we may then perform a projection to derive the single-particle partition functions 
on S'^. Such projection amounts in a consistent truncation of A/" = 4 super Yang Mills as 
discussed in [SU], and it can be realized by taking into account that the only modes that 
actually contribute to the partition function on S*^ are those for which the eigenvalue of 
J3 is equal to half the monopole charge. The projection onto S"^ can thus be performed 
introducing into the A/" = 4 partition functions a f/(l) projection operator of the form 



where 2q is the integer monopole charge of the BPS vacua on S"^ and as a notation we 
shall assume g > 0. 

The projection from S'^ to S"^ rescales the radius of the sphere by 1/2 thus giving an 
52 of radius i? = 1/2. 

C.l Scalars 

Scalars on can be expanded in scalar spherical harmonics Sj^m,m{^) where m and fh 
take the values — j/2, — j/2 + 1, . . . , j/2 — 1, j/2. The energy of a scalar on with 
radius Rsi = 1, conformally coupled to curvature, is E = j + 1. The partition function 
for a scalar on then is 




(C.l) 



E 




(C.2) 



00 j/2 j/2 



00 



x(l + x) 



.seal. 



X- 



.i+i _ 



2^j+l _ 



(C.3) 



j=0 m=—j/2 m=—j/2 



j=0 



where the lower index on z denotes the spacetime dimension. 
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Inserting the projector flC.2p we easily get the partition function for a scalar on S"^. 
The scalar partition function in the presence of a monopole of charge q becomes 



oo j/2 j/2 
j=0 m=-j/2 m.=-j/2 



" de 



(C.4) 



Performing the sums we end up with an integral 

d9 x{l - x^) cos{2q9) 



vr (1 + — 2a; cos i 



that can be easily done and gives 



2q 



_{l-x 



2\2 



(C.5) 



(C.6) 



We can now reintroduce the appropriate dependence on the radius R = Keeping 
into account that the partition function (1C.6P is defined on an with radius R= 1/2, 
to get the one with a generic radius R = 1/fi amounts in simply replacing 



x^ ^x = e"^'^ 



(C.7) 



without having to compute a single determinant. 



C.2 Vectors 

Vectors on can be expanded in vector spherical harmonics V^>fim(^) which belong to 
the representations (ji, J2) = ^) and (ji, J2) = ('^, '^), respectively. The energy 
for both the representations is given by -E = j ' + 1. The partition function on for the 
+ vector component is then 



(i+l)/2 {i-l)/2 



vect 

yx) = 

j=l m=-ij+l)/2m=-(j-l)/2 j=l 



for the — vector component we obviously have the same result z'!^'^^'{x) = z^'^''^'{x) and 
the sum of these two quantities gives the partition function for a vector on 

zrHx) = zif-{x) + zr'-{x) = 

Inserting now the projector (1C.2I) into (1C.8I) and into the analogous one for V~ we get 
for the + and — vector components respectively 

00 (j+l)/2 0-l)/2 



j=l m=-{j+l)/2 rh=-{j-l)/2 
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d9{3 + x^ - 4x cos 9) cos 2q9 

TT fl + — 2xCOS^)2 
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(1 -X2)2 
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1 -x2 



(C.9) 



and 



oo (i-l)/2 (j+l)/2 



j=l m=_(j_l)/2m=-(j+l)/2 

^ (1 + - 4x COS ^ + 2 COS 29) cos 2g^ 



71" ^fl 
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27r 
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For g = this integral gives 



z^^^^'\x, q = 0) 



x2(l +a;2) 



(1 -x2)2 



and for g 7^ 



z'^'''{x, q) = X 



2g 



X^(l + X^ 



2q 
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(C.ll) 
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:i-x2)2 1-X\ 

The hmit g — > is discontinuous, in complete agreement with the computations done in 
appendix IB. 31 . Therefore the sums of the + and — vector partition functions for q ^ 
give 

r 4r2 

^---(x, q) = zl'^ix, q) + zl''-{x, q) = x'"' 



whereas for g = 



(1 -X2)2 



(x, g = 0) = z-",^^- (x, 0) + z",^^- (x, 0) 



l + 2g 



l+x' 



1 - x" 



(C.13) 



4x2 



(1 -X2)2 

Again, with the substitution (1C.7P we immediately get back the results fl4.13|4.14p . 



(C.14) 



C.3 Fermions 



Fermions on S can be expanded in spinor spherical harmonics F^^^iVL) which belong 
to the representations (ji, J2) = (|, and (ji, ^2) = |), respectively. The energy 
for both the representations is given by = j + 1/2. Therefore on 5*^ we get 

c» i/2 0-l)/2 a. 3/2 

i=l m=-i/2?n=-(j-l)/2 i=0 ^ ' 

for the + fermion component and the same result for the — fermion component. The sum 
of these two quantities gives the partition function for a fermion on 



spin/ \ spin./ \ I spin. I 

Z/ X = Z/i X + Z/_ X 



4x3/2 



4 v-^; - ^4+ v-^; ^ ^4- ^ ~ d _ a;)3- 
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Inserting the projector into flC.15|) and into the analogous one for F one gets the 
partition functions for a + or — spinor on S"^ 



spin / 



oo j/2 (i-l)/2 
j=l m=-i/2m=-(i~l)/2 

de 2x^/2 (1 - X COS 6) COS 2qe 



TX 



?2 — 2xcos^)2 



" de 



g2ie(m-g)^j+l/2 



2^3/2 2gx3/2 



_(1-X2)2 1 



(C.16) 



Z_ (x 



oo 0-l)/2 i/2 



(i+1) / ]^ 2ie{fh-q) j + l/2 

jr-=l m=-(j-l)/2 m=-j/2 

^ de 2x^/2 ^ COS e) COS 2g^ _ _^2g 



TT 



2x COS ( 



— ( 

27r 

2x^/2 



X 



2^2 



+ 



gx 



1/2 



X^ 



(C.17) 



Adding (1C.16P and (1C.17P we get the partition function for a fermion on 5*2 in the non- 
trivial background 



^^^-•(a;, q) = "(x, q) + zT'^Xx, q) = x 



2q 



2x2(xi/2 + x"i/2) /xi/2(l + a;) 
~1~ Zq 



X 



2\2 



X^ 



(C.18) 

With the substitution flClTl) we get back the resuh fliASll . 

The complete partition function for our theory can now be constructed using (14.51) . As 
we showed before the presence of the monopole background (13.41) breaks the original U (N) 
invariance to the subgroup 117=1 U{Nj) so that the positive definite charge 2q, appearing 
in the single-particle partition functions, is actually a function of the integers labelling 
the sectors into which the monopole field splits. It can be written here as 



\nj - nj\ 



We easily get 



/k f k oo ^ 

[l[dUj]exp\ E-K'(^") + (-l)"^'4'(^")]Tr(f/7)Tr(([/. 
7=1 l/,J=ln=l 

(C.19) 

Here k is the number of sectors into which the monopole field splits and reintroducing the 
appropriate dependence on the radius R = l/n with the substitution (lC.7p . we recover 
for the bosonic partition function 



zi\x,q) = 6x'?"+^/2 



:i+x) 



2q 



ij 



;i-x) 

and for the fermionic one 

z^p^{x,q) = 4x^" 



X 



+ x' 



qij 



4x 



X 



l + 2q_ 



IJ 



X 



1 — X 



2x(xi/^ + x-^/^) ^ /x^/'^(l + v^) 



(l-x)2 
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X 



(C.20) 



(c.2i; 



We thus reobtain with a very simple and straightforward procedure the result fl4.33p . 
up to the constant (temperature-independent) Casimir contribution. Of course the path 
integral approach has the advantages of giving to Tr(t//) the meaning of matrix holonomy 
along the thermal circle and of providing an explicit derivation of the Casimir energies. 



D Solving the matrix model 

The solution of the matrix model in the presence of a logarithmic interaction has been 
reduced, in section 17. to solve the non-linear differential equation fl7.19p with a given 
set of boundary conditions. Surprisingly, this equation can be explicitly integrated. To 
achieve this goal, we first express p{t) in terms of t and p'{t) by means of (17.191) 

p^'^- iQp'it)mt)-i) ■ ^ ^ ^ 

Subsequently we take the derivative of with respect to t on both sides. The differential 

equation (17.191) factorizes into two factors, which can be set separately to zero. In fact, 
we obtain 

{256tp{tY - I28tp{tf + 16 + t - l) p'{tf - 8pV(t) + p"{t) = 0, (D.23) 
which implies 

p"(t) = ^ p{t) = At + B (D.24) 

and 

256tp'(t)^ - I2%tp'{tf + 16 + t - 1) p'{tf - 8p^p'{t) + = 0. (D.25) 

Consider first (]D.24p . This solution can only satisfy the boundary condition (s) associated 
to the strong-coupling region (see section 17.11) and thus it seems the natural candidate to 
generate T^itjp). This implies that the integration constant B is fixed to be — |p(p + 1). 
The constant A is instead determined by imposing that ( ID. 241) actually solves (I7.19IF ^I. 
We obtain 

vt 1 

The free energy TQ(t,p) is evaluated by integrating (I7.16P with the boundary condition 



•^0 P) = -l ((log(4) - 3)P + (P + 1)' log(p + 1) - log(p)) + - I log(t) . 

(D.27) 

As discussed in section 7.1, (ID. 271) is not the right solution at small t because cannot 
reproduce the series obtained from the large p expansions. We come now to (lD.25p . It 
is an algebraic quartic equation, which determines p'{t) as a function of t and p. We 



"'^■'Since we have taken a derivative of I7.19[ we could have potentially added spurious solutions 
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have four solutions, whose quahtative behavior can be investigated by writing the inverse 
function 

^ ^ (4(p-lMt)-p) [A{l+p)p'{t)-p) 

16p'(t)2 (4p'(t) - 1)' ^ ■ ' 

and by drawing its plot. Since we are interested in positive t and in real solutions, we 
can focus our attention just on the interval [-^^^, The plot is given in fig. [61 We 



4 (p+ 1) 



P 

4 (p-1) 



Figure 6: Plot of the r.h.s of (jD.28p . It diverges for p' = 1/4. For any positive t we have two solutions. 

immediately recognize that there are two potential solutions in this region. At small t, 
they are both finite and their values at t = are respectively 



Pi 



P 



4(p 



and P2 



P 



(D.29) 



For large t, both solutions approach 1/4 but with opposite subleading term. In fact, by 
setting p ~ 1/4 + in flD.281) . we immediately find 



1 
4 



0{t) and p'^it) 



1 

4 



1 



4Vt 



0{t). 



(D.30) 



The actual functions pi,2(^) can be easily recovered by exploiting (ID. 221) . which provides p 
in terms of p' and t (and p). It is easy to check that the solution p2{t) can be dropped since 
its behavior at small and large t is in contrast with the boundary conditions. Instead, 
we can identify pi(t) with the weak-coupling solution pwif) and by integrating ( 17.20p to 
evaluate ^^{t^p). The integration constant is fixed by requiring that our free energy 
coincides with that of the Gross- Witten model for large t. The logarithmic interaction 
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is in fact sub-leading for t ^ 1. Nicely the integration over t can be performed without 
an explicit knowledge of Pw{t)- In fact flD.28|) defines an invertible mapping in the range 
■^^^ < p' < I (see fig. [6]). Thus, by means of flD.281) . we can write 



Pw{t) 



t 

, ^ , ,Ap' i^P'^ - 1)' - 64 ip'J) (/ i^p'^ - if + 16 ip'J (Apl + 1)) 
fw+ / dp^ 



32 (1 - Ap'J ip'j' (p2 (1 - AplY - 16 ip'j') 
/. + ^ (f - ^ + 16 (log (p'J - 2.tanh- (. + 4 (1 - .) p'^) + 



+ log (1 - 4p'J + : ^ 



4PL 



(D.31) 



Here is the arbitrary constant of integration. Requiring that we reobtain the usual 
Gross- Witten model for t ^ 1 fixes our constant to be 

U = -l + \pii-^ + log(16))p - 2 log(p - 1) + 2 log(p + 1)). (D.32) 
With this choice expansion of the free energy J-'Q'{t,p) for large t takes the form 

.2 1 /l\3/2 



^Vl H- i (log (i) - 3) - Ip^yi - 1! . 1-/ (/ - 4) (i)^ 



(D.33) 



The leading behavior is independent of p and it coincides with that of the Gross- Witten 
model. The above expression contains also the result of the semiclassical approximation 
(17.151) . up to higher orders in p'^"' can also compute the small t behavior of 

this solution and it is given by 

^oit^P) =\ (log(p)/ + (3 - log4)p - {p + l)'log(p + 1)) + 

Surprisingly, we see that J^Q(t,p) satisfies also the boundary condition (17. 8p for small t 
and reproduces, in that regime, the result of the large p expansion. In other words, flD.311) 
and flD.321) provide a solution which smoothly interpolates between the strong and the 
weak coupling regime. 
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